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1.0  Abstract 

Long-term  monitoring  (LTM)  costs  make  up  a  substantial  percentage  of  environmental 
restoration  budgets  for  both  government  agencies  and  the  private  sector.  A  major  portion  of 
LTM  costs  are  costs  associated  with  groundwater  sampling  and  analysis.  One  promising  cost¬ 
saving  strategy  is  the  replacement  of  traditional  groundwater  sampling  approaches  by  a 
combination  of  passive  sampling  devices  and  in  situ  sensors.  However,  regulatory  and 
operational  acceptance  of  changes  in  groundwater  sampling  methodologies  is  conditional  on 
both  laboratory  and  field  validation  of  the  methods  and  an  in-depth  understanding  of  the  way 
such  methods  represent  groundwater  chemistry. 

Objective.  The  project  objective  for  ER-1704  was  to  develop  a  fuller  understanding  and 
description  of  how  contaminant  concentrations  measured  in  a  well — ^using  either  passive 
sampling  devices  or  in  situ  sensors — relate  to  contaminant  concentrations  in  the  surrounding 
formation.  ER-1704  research  has  elucidated  several  key  factors  in  the  relationship  between 
solute  concentrations  in  wells  and  in  the  surrounding  formation,  thus  providing  additional 
scientific  basis  for  cost  effective  application  of  alternative  sampling  and  monitoring  strategies. 
Solid  scientific  basis  for  non-purge  groundwater  monitoring  is  imperative  for  gaining  user  and 
regulatory  acceptance  of  the  alternate  techniques.  Defining  the  measurable  attributes  of 
aquifer/well  conditions  that  either  suggest  likelihood  or  imply  limits  to  well  mixing  phenomena 
were  primary  objectives  of  the  work. 

Methods.  Physical  and  numerical  modeling  efforts  were  completed  to  illustrate  the  relationship 
between  solutes  in  groundwater  and  those  found  in  wells.  Physical  mixing  processes  that  occur 
during  ambient  flow  in  aquifers  and  during  active  pumping  activities  were  explored  to  show  how 
potentially  stratified  contaminants  may  manifest  in  open  well  screens.  Complex  flow  patterns 
driven  by  flow  velocity  changes,  very  small  density  contrasts,  and  temperature  contrasts  were 
identified  as  possible  sources  of  well  convection  and  mixing.  Specifically,  the  modeling 
explored; 

•  2D  and  3D  physical  tank  models.  Dye  tracer  testing  was  conducted  to  illustrate  flow  in 
and  around  models  wells. 

•  Analytical  models.  Models  of  well  flow  patterns  were  used  to  resolve  well  flushing 
during  purging. 

•  Axisymmetric  Numerical  Models.  Pumping  models  were  constructed  to  simulate 
pumping  flow  to  a  well,  allowing  variations  in  adjacent  hydraulic  properties  and 
contaminant  distribution  to  be  easily  tested. 

•  Numerical  Bore  Model.  In-well  density  convection  phenomena  were  explored  using  a 
closed  bore  well  model. 
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•  Numerical  Section  Model.  2D  mass  and  energy  transport  models  were  constructed  to  test 
inflow  phenomena,  density-driven  convection,  and  the  effects  of  in-well  flow  limiting 
devices. 

Physical  field  testing  was  completed  to  test  and  illustrate  contaminant  fluctuations  and 
differentials  in  field  conditions.  As  such,  these  tests  were  less  well-controlled  and  near-well  flow 
conditions  were  only  inferentially  known  until  experimental  data  was  generated.  This  approach 
is  much  more  similar  to  practitioner  conditions  where  collected  data  drives  interpretation. 

•  Multilevel  passive  sampling.  Multilevel  sampling  was  conducted  using  two  approaches; 
1)  open  well  sampling;  2)  isolated  zone  sampling. 

•  Well  Installation.  New  wells  were  installed  with  special  design  configurations  to  explore 
the  relationship  between  contaminant  concentrations  inside  and  outside  the  well  during 
ambient  flow  and  during  pumping. 

•  Purge  Dynamics  Testing.  Special  purging  protocols  tested  how  pumping  rate, 
consistency  of  flow  rate,  changes  in  flow  rate  prior  to  sampling,  and  duration  of  pumping 
affects  pump  discharge  concentrations,  relative  to  concentrations  measured  in  the  well. 

•  Bottle  Filling  Test.  An  outdoor  field  test  of  methods  that  are  commonly  used  to  fill  VGA 
vials  was  conducted.  This  test  was  conducted  to  assess  the  degree  to  which  field  sample 
handling  affects  analytical  results. 

Findings. 

1.  Passive  sampling  (or  remote  passive  sensing)  in  unpumped  wells.  Wells  tend  to 
redistribute  adjacent  aquifer  contaminant  heterogeneity.  Homogenization  often  occurs, 
but  the  degree  of  mixing  varies  from  well  to  well.  Redistributive  effects  are  highly 
sensitive  to  very  small  density  contrasts  in  the  inflowing  fluid.  Because  of  these  density- 
driven  effects,  when  aquifer  contaminant  chemistry  is  stratified,  solute  distribution  in  a 
well  may  not  match  that  in  the  adjacent  formation.  Often  times,  the  passive  sample 
closely  matches  the  flow-weighted-average  of  inflowing  water.  This  condition  is  not 
perfect,  however,  and  concentration  distributions  vary  from  well  to  well.  There  are  many 
factors  affecting  contaminant  distribution  in  an  unpumped  well  that  are  not  readily  known 
without  exhaustive  testing.  Pumping  dynamics  are  also  complex,  meaning  that 
differences  from  pumped  comparators  do  not  always  clarify  which  method  is  “correct”  in 
such  circumstances.  Contaminant  position  in  the  aquifer  relative  to  the  passive  sampling 
device  position  in  the  well  is  usually  not  known,  and  understanding  of  that  aquifer- 
contaminant  geometry  may  not  be  improved  by  multilevel  passive  sampling  in  an  open 
bore.  Isolated  zone  sampling  can  improve  determination  of  aquifer  contaminant 
stratification,  but  the  degree  of  improvement  is  variable  and  well-dependent.  As  a 
general  replacement  for  purge  sampling,  passive  sampling  or  passive  sensing  do  yield 
similar  results  in  most  cases.  Where  results  are  different,  and  information  is  desired. 
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causes  can  often  be  found  through  more  thorough  multiple-sample-testing  during 
purging,  multilevel  passive  sampling  within  the  well,  or  a  combination  of  the  two. 

2.  Pumping  dynamics.  Aquifer  eontaminant  position  relative  to  the  pump  intake  position  is 
rarely  known,  but  knowledge  may  be  improved  by  colleeting  multiple  samples  during 
purging.  Contaminant  position  in  the  aquifer  relative  to  the  pump  position  in  the  well 
drives  concentration  stability  during  pumping.  Volume  of  water  removed  from  the  well 
is  a  more  reliable  predictor  of  contaminant  concentration  stability  than  measurement  of 
traditional  “purge  parameters”  such  as  temperature  pH,  EC,  ORP  and  DO;  yet 
concurrently,  large  volumes  removed  do  not  always  assure  contaminant  concentration 
stability.  Some  contaminants  of  interest  stabilize  faster  than  others,  depending  on  the 
specific  well,  suggesting  that  unique  contaminant  distribution  heterogeneity  and 
chemical-specific  biological  activity  can  influence  the  stability  of  contaminant 
concentrations  during  purging. 

3.  Physical,  Numerical,  and  Analytic  Models.  Good  matches  between  simulation  results, 
experimental  data,  and  theoretical  analysis  of  flow  and  transport  support  our  hypothesis 
that  redistributive  effects  by  vertical  mixing  is  common.  Physical  model  experiments 
show  that  small  heterogeneities  exert  strong  influence  on  flow  in  the  open  bore  of  a 
model  well.  Horizontal  laminar  flow  across  the  model  bore  could  not  be  reproduced 
except  under  conditions  in  which  density  variations  are  much  less  than  would  be 
expected  in  the  field.  Density  contrasts  equivalent  to  as  little  as  10  parts  per  billion  of 
dissolved  solids  are  enough  to  cause  near  complete  vertical  redistribution  in  simulated 
wells.  Thermal  convective  behavior  indueed  by  shallow  seasonal  thermal  changes  and 
deeper  geothermal  gradients  may  also  either  promote  or  inhibit  well  mixing  effects. 

Overall  Conclusions  and  Technical  Transfer.  Contaminant  redistributive  effects  in  wells  are 
nearly  always  present.  Complete  mixing  appears  to  be  very  common;  however,  it  is  not 
universal.  There  is  a  continual  balance  between  inflowing  contaminant  stratification  (where 
present)  and  factors  driving  in-well  mixing.  Findings  here  imply  common  and  very  small  drivers 
are  responsible  for  slow  but  vigorous  mixing  relative  to  the  residence  time  of  water  flowing 
through  a  typieal  well  screen.  Therefore,  a  tendency  toward  homogenization  is  anticipated  to  be 
common  in  field  eonditions.  Most  wells  should  experienee  strong  redistribution  effects,  but 
some  wells  may  maintain  stratification  or  perhaps  re-stratify  differently  from  the  surrounding 
formation.  Ongoing  technical  transfer  of  these  findings  will  promote  better  understanding  in  the 
environmental  community  that  wells  often  represent  a  mixed  flow-weighted  average  of  the 
adjacent  formation  chemistry.  This  better  understanding  will  yield  cost  savings  in  both  short¬ 
term  and  long-term  timeframes  by  accelerating  the  approval  process  for  non-purge  alternative 
sampling  strategies,  including  passive  sampling  and  in  situ  sensor  technologies. 


XV 
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2.0  Objective 

Regulatory  acceptance  of  in  situ  groundwater  monitoring  tools,  as  replacements  for  traditional 
groundwater  sampling  methods,  is  conditional  on  confidence  in  and  an  in-depth  understanding  of 
how  such  tools  provide  accurate  measurements  of  groundwater  chemistry.  Providing  this 
confidence  and  understanding  was  a  principal  objective  of  SERDP  Statement  of  Need  (SON)  09- 
05.  The  ER-1704  project  plan  incorporated  efforts  designed  to  explain  the  mechanics  of  ambient 
groundwater  flow  in  and  around  monitoring  wells  and  during  active  sampling.  Efforts  during  the 
project  have  yielded  promising  results  that  help  explain  these  mechanics.  Deeper  understanding 
and  technical  transfer  of  that  understanding  is  essential  for  regulatory  acceptance  of  these 
alternative  approaches  to  long-term  monitoring  (ETM)  for  groundwater. 

Considerable  inertia  still  exists  within  the  groundwater  sampling  community  with  regard  to  a 
perceived  need  to  purge  groundwater  wells  prior  to  sampling  them.  While  this  inertia  has  been 
overcome  to  a  degree  by  efforts  of  the  project  team  members  in  their  work  with  the  Interstate 
Technical  and  Regulatory  Council  (ITRC)  Passive  Sampler  Team,  American  Society  for  Testing 
and  Methods  (ASTM),  and  other  groups,  there  is  substantial  work  to  be  done  to  allow  passive 
methods  to  be  widely  and  fully  accepted  where  appropriate.  Since  in  situ  sensor  technology  is 
also  by  definition  passive,  application  of  in  situ  sensors  also  implicitly  requires  acceptance  of 
passive  sampling  techniques.  Our  approach  focused  on  the  enhanced  understanding  of 
representativeness  of  passive  sampling  methodologies.  This  approach  is  designed  to  provide  a 
basis  for  avoiding  the  “purge,  sample,  analyze”  standard  protocol  that  is  largely  followed  today. 

Accordingly,  ER-1704  project  work  was  designed  to  develop  a  fuller  understanding  and 
description  of  how  contaminant  concentrations  measured  in  a  well  relate  to  contaminant 
concentrations  in  the  surrounding  formation.  Specifically,  the  project  addresses  the  following 
questions; 

1.  How  does  the  relationship  between  solute  concentrations  in  the  well  and  in  the 
surrounding  formation  depend  on  geochemical  conditions,  subsurface  heterogeneity,  well 
construction,  sampling  methodology  and  other  primary  controls  on  solute  transport  in  the 
subsurface? 

2.  Several  lines  of  evidence  now  indicate  that  rapid  vertical  mixing  occurs  in  wells  even 
under  ambient  (non-pumping)  flow  conditions.  Which  mechanism(s)  cause  this 
enhanced  mixing? 

3.  Can  results  from  previous  studies,  combined  with  experimental  studies  and  modeling 
studies  performed  for  the  proposed  project,  be  used  to  develop  a  comprehensive 
conceptual  model  of  well  behavior,  and  thus  be  used  as  a  guide  to  developing  a  reliable 
predictive  model  for  new  subsurface  monitoring  methods? 

4.  Based  on  an  improved  understanding  of  well  behavior,  what  adaptations  might  be  made 
to  existing  sampling  methods  to  provide  more  and/or  more  accurate  information  about 
subsurface  conditions? 
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To  meet  the  technical  objectives  of  this  project,  we  proposed  a  technical  approach  that  included 
the  following  components. 

(1)  Formulation  of  a  comprehensive  conceptual  model  for  flow  and  transport  into  and  around 
a  well  which  incorporates  all  relevant  physical  properties  and  processes. 

(2)  Laboratory  experiments  in  2D  and  3D  sand-tank  physical  models,  focused  on  providing 
experimental  data  allowing  for  numerical  model  development  and  calibration.  In  this 
objective  we  expanded  on  the  work  done  by  (Britt  2005)  through  a  combination  of  2D 
and  3D  experiments  that: 

o  specifically  address  the  mechanisms  that  lead  to  rapid  vertical  mixing  in  a  well, 

o  explore  how  a  range  of  heterogeneities,  both  physical  and  geochemical,  affect 
concentrations  in  the  well, 

o  explore  how  different  pumping  rates,  and  changes  in  pumping  rate,  affect 
concentration  distribution  in  the  well, 

o  examine  how  non-equilibrium  chemistry,  combined  with  other  transient  effects 
(e.g.  pumping  effects)  relate  to  concentrations  measured  in  the  well,  and 

o  expand  previous  physical  model  studies,  which  essentially  simulated  a  2D  flow 
problem,  to  examine  3D  behavior  around  a  well,  in  order  to  extend  conclusions 
based  on  the  2D  studies  to  actual  well  behavior. 

(3)  Application  of  numerical  models  of  flow  and  solute  transport  in  and  around  a  well  to 
simulate  the  experiments  conducted  in  component  2.  These  models  are  based  on  the 
conceptual  model  developed  in  component  1.  The  numerical  models  will  enhance  our 
understanding  of  the  experimentally  observed  phenomena,  as  well  as  establish  and 
sharpen  our  ability  to  accurately  predict  flow  and  transport  into  and  around  a  well.  This 
work  builds  on  previous  results  of  physical  and  numerical  modeling  by  other  scientists 
(e.g.  Martin-Hayden  2000a,  Varljen  et  al.  2006)  to  arrive  at  a  comprehensive  predictive 
model  which  accounts  for  the  effects  of  the  primary  controls  on  contaminant 
concentration.  The  focus  of  this  component  will  be  on  incorporating,  based  on  a  sound 
conceptual  model,  the  effects  of  in-well  behavior  (stratification  and  mixing)  and 
geochemical  processes  in  these  previously  developed  3D  models. 

(4)  Field  experiments  using  different  sampling  methods  focuses  on  field  testing  and 
validation  of  the  predictive  capabilities  of  the  numerical  models.  The  field  work  includes 
experiments  using  both  passive  and  active  sampling  methods  to  evaluate  the  field-based 
predictive  capability  of  our  models.  Multiple  downhole  sampling  approaches  are  utilized 
to  isolate  differences  found  between  methods.  These  methods  include  open-well  interval 
sampling,  after  Parsons  (2003);  baffle  isolation  zone  testing,  after  Britt  and  Calabria 
(2008);  purge  dynamics  testing,  after  Martin-Flayden  and  Britt  (2006);  and  other  methods 
described  below. 
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3.0  Background 

A  brief  overview  of  well  sampling  approaches  with  respect  to  well  flow  dynamics  is  provided 
here  (adapted  and  updated  from  Britt  (2005)).  Some  of  the  unknowns  identified  therein  were 
investigated  in  the  (2005)  study,  and  in  the  present  investigation  in  response  to  the  SERDP  SON. 

Purging  monitoring  wells  in  order  to  collect  representative  samples  has  been 
conventional  procedure  for  upwards  of  25  years  (e.g.  USEPA,  1977,  EISGS,  1980).  Over 
much  of  that  period,  some  level  of  controversy  has  continued  over  how  best  to  purge 
wells,  or  whether  it  is  necessary  to  purge  wells  at  all  (Robin  and  Gillham,  1987;  Powell 
and  Puls,  1993,  Newell  et  ah,  2000).  Many  studies  have  identified  how  anomalous  or 
otherwise  unrepresentative  results  may  be  generated  from  traditional  purge  and  sampling 
techniques  (e.g.  Robbins,  1989,  Gibs  and  Imbrigiotta,  1990;  Reilly  and  Gibs,  1993; 
Conant  et  ah,  1995;  Church  and  Granato,  1996,  Martin-Hayden  and  Robbins,  1997, 
Reilly  and  EeBlanc,  1998,  Hutchins  and  Acree,  2000,  Elci,  et  ah,  2001  and  2003). 
Common  situations  shown  by  these  and  other  studies  have  caused  investigators  to 
question  what  traditional  purge-and-sample  ground  water  monitoring  results  represent. 
Low-flow  purge  techniques  (e.g.  Puls  and  Barcelona,  1996)  were  developed  and  widely 
adopted  to  address  some  of  the  problems  (e.g.  falsely  elevated  metals  concentrations), 
with  the  added  benefit  of  reducing  purge-water  waste.  No-purge  techniques  are  also 
being  explored  and  are  being  adopted  where  applicable  (Vroblesky,  2001,  Parker  and 
Clark,  2002;  Parker  and  Mulherin,  2007,  Parsons  Engineering  Science,  2003;  ITRC, 
2004,  2006,  2007).  These  alternative  techniques  can  solve  some  problems  like  elevated 
turbidity  and  VOC  loss  caused  by  bailer  agitation  or  high  pump  rates,  but  they  do  not 
solve  problems  with  vertical  flow  (Elci,  et  ah,  2001  and  2003),  pumping-induced 
variability  (Martin-Hayden,  2000a  and  2000b;  Gibs,  et  ah,  2000,  Martin-Hayden  and 
Britt  2006),  or  well  convection  (Martin-Hayden  and  Britt  2006;  Vroblesky  et  al.  2007). 
Other  techniques  involve  installing  multichannel  tubing  wells  (Einarson  and  Cherry, 
2002);  short-screen  direct  push  wells  (Kram,  et  ah,  2001);  or  they  involve  installation  of 
devices  such  as  the  Discrete  Multi-Level  Sampler  within  existing  longer  screen  wells 
(Puls  and  Paul,  1997),  or  installation  of  baffle-type  devices  (Vroblesky  el  al.  2007,  Britt 
2008).  These  techniques  are  effective  but  may  require  a  long  term  commitment  to 
relatively  expensive  multi-interval  sampling. 

Contaminant  transport  interpretation  and  cleanup  decision-making  directly  relies  on 
accurate  data  and  a  clear  understanding  of  what  the  data  represent.  Misdirection  from 
anomalous  data  should  clearly  be  avoided.  Church  and  Granato  (1996),  Reilly,  et  al. 
(1989),  and  Elci,  et  al.  (2001)  show  the  mechanics  of  sampling  bias  and  misdirection 
caused  by  vertical  well  bore  flow.  However,  current  work  by  this  author  and  colleagues 
includes  feasibility  assessment  of  long-screen  well  rehabilitation  by  physically  limiting 
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vertical  flow  within  wellbores  as  a  partial  solution  to  vertical  ambient  flow  (Britt  and 
Calabria  2008). 

Other  sources  of  misdirection  from  sampling  data  include  partial  penetration  (Robbins, 
1989;  Martin-Hayden  and  Robbins,  1997;  and  Chaing  et  al.,  1995)  and  lithologic 
heterogeneity  (Gibs  et  al.,  1993).  Anderson  and  others  (1992)  showed  how  DNAPL 
fingers  and  pools  result  in  very  strong  dissolved  contaminant  stratification  in  an  aquifer. 
Well  screen  dilution  and  redistribution  of  strong  contamination  stratification  could 
clearly  misdirect  an  investigator  attempting  to  delineate  contaminant  distribution  and 
transport.  But  very  short  multi-interval  wells  can  also  miss  vertically  limited,  highly 
contaminated,  zones.  Just  as  the  environmental  industry  is  moving  toward  continuous 
geologic  logging  and  chemical  mapping  of  boreholes  (e.g.  GeoVis,  MIP  [membrane 
interface  probe])  rather  than  collecting  soil  samples  once  every  five  feet,  ground  water 
data  should  similarly  be  collected  over  entire  intervals.  Of  course,  this  is  limited  by 
geology,  and  wells  should  not  cross  significant  lithologic  boundaries,  but  if  one  doesn’t 
know  where  contaminants  are  moving  within  a  defined  aquifer  zone,  one  should  sample 
and  characterize  the  entire  interval.  Ground  water  monitoring  system  design  must 
incorporate  solutions  for  these  conflicting  problems  to  generate  the  most  useful  sampling 
data. 

To  alleviate  concerns  about  these  conflicting  problems  with  representativeness,  multi¬ 
level  sampling  within  existing  wells  (Puls  and  Paul,  1997;  Vroblesky  and  Peters,  2000; 
Parsons  Engineering  Science,  2003;  Britt,  2006;  Britt,  2008;  Britt  and  Calabria,  2008) 
and  sampling  multi-interval  well  installations  (E inarson  and  Cherry,  2002,  and  Kram,  et 
ah,  2001)  can  certainly  help,  as  long  as  sufficient  effort  is  made  to  understand  what  the 
results  indicate,  and  one  can  afford  to  generate  consistent  comparison  data. 

Notwithstanding  the  interest  of  multi  level  sampling,  the  practical  lons-term  objective 
for  most  sampling  programs  is  to  be  able  to  efficiently  collect  a  single  representative 

sample  from  each  well  The  working  hypothesis  for  the  (2005)  study  was  that  under 
simple  scenarios,  contaminant  mixing  and  dilution  can  and  does  occur  in  monitoring 
wells  under  natural  flow-through  conditions.  The  flow-through  aspect  of  this  hypothesis 
is  no  longer  considered  contentious — as  ASTM  observes  in  its  2002  standard  practice 
document  for  low-flow  purging; 

“Low-flow  purging... is  based  on  the  observations  of  many  researchers 
that  water  moving  through  the  formation  also  moves  through  the  well 
screen.  Thus,  the  water  in  the  well  screen  is  representative  of  the 
formation  surrounding  the  screen.”  (ASTM,  2002) 

If  flow-through  and  mixing  creates  a  flow-weighted  average  concentration  in  the  screen 
interval  of  a  well,  it  is  arguable  that  wells  are  “naturally  purged.”  Simple  tracer 


ER-1704  Final  Report  2014 


dissipation  tests  (see  Robin  and  Gillham,  1987)  could  determine  whether  individual  wells 
are  “naturally  purged”  or  not.  Consequently,  naturally  pursed  wells  could  allow 
collection  of  the  ideal  sample — a  sinsle,  inexpensive,  representative  sample  collected 

directly  from  the  screen  interval  of  the  well. 
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4.0  Materials  and  Methods 

This  project  was  designed  to  evaluate  whether  naturally  purged  wells  can  proved  an  “ideal 
sample,”  and  to  explore  in  more  detail  the  dynamics  of  how  measured  contaminant 
concentrations  vary  as  a  function  of  parameters  such  as  hydrogeologic  conditions, 
geochemistry,  well  type,  sampling  method,  contaminant  type  and  concentration.  Physical 
experiments  included  to  quasi  2-dimensional  models  of  near  well  environments.  Both  porous 
media  and  simulated  fracture  media  were  physically  modeled.  Numerical  models  are 
explained  in  greater  detail  as  they  are  more  innovative  and  unique  from  a  materials  and 
methods  perspective.  The  physical  models  are  somewhat  intuitive  therefore  not  explained  in 
great  detail  here  beyond  dimensions  of  the  experiments.  Britt  (2005)  and  Martin  Hayden  and 
Wolfe  (2000)  described  similar  physical  “sand  tank”  models  in  more  detail.  Methods  for  this 
approach  follow  below. 

4.1  Conceptual  Model  Development 

Existing  conceptual  models  for  flow  and  transport  around  wells  assume  Darcian  flow  to  and 
through  the  well,  which,  in  turn,  implies  laminar  flow  to  and  through  the  well  screen.  However, 
prior  to  the  experiments  in  this  project,  field  data  (Britt  and  Calabria  2008)  and  experimental  data 
(Britt  2005)  suggested  that  laminar  flow  through  the  well  does  not  hold  in  tested  laboratory 
conditions  and  limited  testing  in  field  conditions.  In  this  component  modifications  and 
enhancement  to  existing  conceptual  models  were  developed  necessary  to  reconcile  observations 
with  theory.  This  aspect  of  the  research  was  an  iterative  effort,  led  by  the  results  from 
experimental  laboratory  tests  using  2D  and  3D  physical  models,  analytical  and  numerical 
modelling,  and  field  tests. 

4.2  Physical  Laboratory  Experiments 

Previous  experiments  (Britt  2005,  Figure  1)  demonstrated  that  under  ambient  conditions,  vertical 
mixing  in  wells  is  much  more  rapid  than  previously  anticipated.  Further,  the  horizontal  flow 
assumption  could  not  be  replicated.  Faboratory  experiments  for  this  project  follow  up  that  work, 
by  examining  the  same  effect  under  different  conditions,  and  determining  how  well  construction 
affects  the  rate  of  vertical  mixing.  In  addition,  the  following  issues  related  to  well  behavior 
under  different  sampling  conditions  in  different  environments  are  explored. 

•  Fractured  media  transport  and  flow  mixing  in  a  model  well 

•  Pumping  flow  dynamics  testing  to  establish  comparison  criteria  to  traditional  sampling 


Figure  1.  Sand  tank  apparatus  (A)  used  in  Britt’s  (2005)  dye  tracer  experiments  in  a  quasi-  2D  horizontal  flow  field  and 
(B)  a  time  lapse  sequence  of  images  illustrating  rapid  vertical  movement  of  dye  relative  to  horizontal  movement.  Dye 
was  injected  in  the  sand  to  the  left  of  the  well  and  entered  the  well  as  illustrated  in  panel  A.  Times  beneath  time  lapse 
images  denote  hours  from  time  that  dye  was  first  observed  in  the  well.  Dye  solution  in  this  experiment  was  very 
slightly  denser  than  the  ambient  flow  solution. 


4.2. 1  2D  Physical  Modeling 

Quasi-two-dimensional  ^  physical  models  were  developed  to  examine  a  range  of  physical 
processes  for  which  the  3D  flow  regime  around  the  well  is  considered  subordinate  to  questions 
about  the  interaction  of  the  formation  with  the  slow  moving  water  in  the  well.  For  these 
experiments,  a  “model  well”  is  defined  as  the  open  portion  of  the  physical  system  than  comprises 
the  overall  physical  model.  This  is  in  contrast  to  the  aquifer  portion  of  the  model,  which  may  be 
fdled  with  porous  media,  or  a  simulated  fracture  system.  The  “model  well”  is  designed  to 
replicate  field  conditions  of  an  open  section  of  a  well  or  borehole  within  an  aquifer. 

•  Effect  of  well  construction  on  internal  mixing 

•  Fractured  media  transport  and  flow  mixing  in  a  model  well 

•  Strong  hydraulic  conductivity  contrast  flow  testing  in  a  model  well 

•  Flow  dynamics  testing  to  establish  comparison  criteria  to  traditional  sampling 

Physical  models  were  constructed  to  allow  rapid  change  out  of  filter  and  well  screen  material  to 
accommodate  different  testing  scenarios.  Tests  were  conducted  with  sufficient  replications  to 
satisfy  the  investigatory  team  that  the  observed  phenomena  are  repeatable  and  reliable  for 
predictive  purposes.  Observations  were  recorded  using  time-lapse  photography,  data  logging, 
3D  planar-laser-induced  fluorescence  tomography,  and  other  methods  depending  on  the  specific 
experiment. 

Individual  test  runs  under  a  variety  of  scenarios  illustrate  an  in-well  mixing  effect  similar  to 
Britt’s  2005  results,  with  additional  details  and  causality  defined  during  the  trials.  Attempts  were 


^  te.  the  model  is  designed  to  produce  a  condition  in  which  flow  and  transport  is  largely  restricted  to  two 
dimensions 
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made  to  find  scenarios  and  conditions  that  did  not  replicate  the  previously  observed  mixing 
phenomena.  The  purpose  of  this  phase  of  the  investigation  is  to  add  to  the  predictive  capability 
and  identification  of  conditions  contrary  to  the  baseline  flow-weighted  averaging  effect,  and 
ultimately  the  project  goal — delineation  of  when  a  single  passive  sample  /  sensor  approach  will 
and  will  not  satisfactorily  provide  a  representation  of  a  well  screen  interval. 

4.2.2  3D  planar-laser-induced  fluorescence  tomography 

We  employed  laser- induced  fluorescence  tomography  to  obtain  a  3D  description  of  the  migration 
of  the  dye  within  the  well  and  to  develop  a  means  of  directly  measuring  dye  concentrations  at 
concentrations  difficult  to  detect  visibly.  Figure  2  and  3  illustrate  the  tank  configuration  and 
apparatus.  A  100  mW  532-nm  wavelength  laser  was  aimed  at  a  cylindrical  lens  to  produce  a 
laser  sheet  which  was  then  collimated  via  a  parabolic  mirror  and  directed  into  the  top  of  the  well 
section  of  the  sand  tank  model.  A  pivoting  dielectric  mirror  was  used  to  translate  the  laser  sheet 
through  the  well  and  images  of  fluoresced  dye  from  a  sequence  of  planes  are  collected  via  digital 
video.  The  laser  plane  translation  rate  was  set  so  as  to  traverse  the  well  column  in  approximately 
one  second,  so  that  each  video  effectively  represents  a  3D  snapshot  of  dye  position  within  the 
well.  The  individual  frames  of  the  1-s  video  scans  were  then  converted  to  3D  data  using  Matlab 
image  processing  functions,  and  the  results  used  to  examine  dye  behavior  during  transport  into 
and  through  the  well. 
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Figure  2.  Photograph  of  INL's  sand  tank,  (A)  before  filling,  and  (B)  a  close-up  of  a  portion  of  the  tank 
during  one  of  the  dye  injection  experiments,  illustrating  aquifer  arrangement  and  dye  movement  in 
the  well  section  of  the  tank 


2014 


ER-1704  Final  Report 


4.2. 3  3D  Physical  Modeling 

Three-dimensional  physical  models  were  developed  to  further  investigate  the  2D  modeling 
findings  in  a  3D  context.  There  are  ongoing  questions  regarding  the  transferability  of  2D 
experimental  work  to  three  dimensions,  and  from  there  to  field  conditions.  Conversely,  the  3D 
models  also  allow  field  findings  to  be  further  explored  under  controlled  laboratory  conditions. 
Testing  the  2D  findings  with  3D  physical  modeling  and  field  experiments  are  critical  to  ultimate 
acceptance  of  the  findings  by  the  environmental  community. 


Figure  3.  Photograph  and  schematic  of  the  planar-laser-induced  fluorescence  tomographic  imaging 
apparatus  developed  to  examine  dye  transport  into  INL's  replication  of  Britt's  2D  sand  tank  model. 


There  are  different  ways  to  approach  3D  physical  model  construction.  Models  can  be  small- 
scale  representations  of  individual  wells  in  cross  section,  or  larger  scale  fully  three  dimensional 
configurations.  There  are  advantages  and  disadvantages  of  each.  A  cross  section  of  a  well 
allows  direct  viewing  and  time-lapse  photography  of  tracer  dyes  during  well-bore  flow/mixing 
and  evaluation  of  downgradient  “shadow”  effects  due  to  that  in-well  mixing,  while  fully  three 
dimensional  models  allow  more  robust  analyses  of  3D  flow  fields  within  and  around  the  well  in 
the  absence  of  transparent  face-plate  effects. 

The  3D  physical  model  was  constructed  to  include  both  the  “Half  Well”  configuration  and  the 
“Full  3D  Well”  in  order  to  submit  both  to  the  same  hydraulic,  solutal  and  thermal  conditions. 
Figure  4  shows  the  two  configurations;  the  Half  Well  preserves  ability  to  visually  track  solute 
movement  along  the  centerline  of  flow  approaching  the  well,  within  the  well  and  down-gradient 
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of  the  well  while  the  Full  3D  Well  allows  visualization  from  the  top  and  bottom  of  the  well 
(through  the  transparent  bottom  of  the  tank). 

4.2.4  3D  Flow  Field 

In  order  to  provide  easy  observations  of  dye 
eoneentration  throughout  their  physieal 
model,  Britt  (2005)  designed  a  sand  tank 
model  that  approximates  a  flow  system 
where  the  well  is  a  continuous  barrier  to 
flow  in  a  plan  view  of  the  system.  Martin- 
Hay  den  and  Wolfe  (2000)  designed  physical 
model  of  a  monitoring  well  cross  section  in 
order  to  track  flow  regimes  and  mixing 
within  the  well  bore.  The  physical  model  in 
this  study  was  designed  to  combine  the 
visualization  of  flow  through  (of  Britt  2005) 
and  quasi-3D  flow  visualization  (of  Martin- 
Hayden  2000b).  The  first  series  of  tests  in 
the  Half  Well  (Figure  4),  represent  a  cross- 
sectional  half  of  a  3-D  flow  field  around  a 
well  in  order  to  directly  observe  the  flow 
regime  and  mixing  within  and  adjacent  to 
the  well  through  the  transparent  face-plate. 

These  tests  were  extended  to  the  Full  3D 
Well  by  using  resistivity  logging  and  in-well  sensors  to  observe  tracer  movement.  These  models 
are  used  to  test  differences  in  behavior  between  the  2D  and  3D  physical  models,  so  that 
conclusions  based  on  the  2D  model  studies  can  be  appropriately  generalized  to  the  3D  case. 

4.2.5  Fracture  Borehole  Model 

Flow  and  mixing  within  fractured  rock  boreholes  are  expected  to  involve  various  processes  that 
distinguish  them  from  monitoring  wells  in  granular  media.  In  order  to  investigate  these 
processes  a  physical  model  of  a  borehole  within  a  fractured  medium  was  constructed  (Figure  5). 
The  borehole  model  consists  of  a  47mm  ID  polycarbonate  tube  intersected  by  two  fractures  each 
consisting  of  a  space  between  two  340mm  x  300mm  plates  of  12.5mm  thick  polycarbonate.  That 
fracture  aperture  is  maintained  with  a  rubber  gasket  around  the  edges  that  can  be  adjusted 
between  0.5mm  and  0.8mm  using  a  set  of  bolts  around  the  perimeter  of  the  plates.  A  hydraulic 
gradient  and  uniform  flow  is  introduced  by  means  of  constant  head  reservoirs  at  either  end  of  the 
fracture.  Dye  is  introduced  at  the  low-head  side  of  the  fracture  and  is  imaged  using  a  DSLR 
camera.  The  lower  fracture,  250mm  from  the  bottom  of  the  borehole,  is  horizontal  and  the  top 
fracture,  570mm  from  the  bottom,  is  inclined  at  a  27°  angle.  The  borehole  extended  another 


T3 


Figure  4.  Physical  model  of  “Half  Well”  and  “Full  3D 
Well”  with  controllable  horizontal  and  vertical  gradients. 
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370mm  above  the  top  fracture  to  represent  the  blank  casing  or  the  unfractured  top  portion  of  the 
borehole. 

Imaging  the  dye  passing  from  the  fracture  into  the  borehole  allows  visualization  of  the  mixing  as 
the  water  enters  the  well  and  quantification  of  the  mixing  within  the  well.  The  zone  down 
gradient  of  the  borehole  represented  the  shadow  zone  indicating  the  influence  of  that  mixing  on 
the  groundwater  in  the  fracture.  The  grid  beneath  the  transparent  fracture  is  used  to  map  the  flow 
regime  as  groundwater  converges  on  the  borehole  (capture  zone)  and  as  well  water  diverges  from 
the  borehole  on  the  down-gradient  side  (shadow  zone).  The  capture  zones  and  shadow  zones  for 
each  fracture  are  expected  to  change  in  geometry  and  size  depending  on  the  head  difference 
between  the  fractures.  The  head  difference  between  the  fractures  is  also  expected  to  control  the 
characteristic  of  mixing  and  transport  within  the  well.  When  density  driven  flow  dominates  the 
flow  within  the  well,  that  mixing  will  determine  the  distribution  of  solutes  in  the  well.  When 
borehole  flow  overcomes  density  mixing  the  borehole  flow  regime  will  dominate  transport. 
Submerging  the  fracture  model  in  a  bath  of  water  maintained  a  uniform  and  constant  temperature 


Figure  5:  Geometry  of  the  fracture  borehole  model,  (a)  downward  oblique  view  of  the  top  angled  fracture,  (b)  a 
side  view  of  the  fracture  model  showing  the  vertical  borehole,  the  angled  top  fracture  and  the  horizontal  bottom 
fracture. 
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as  a  means  of  controlling  thermal  convection  and  density  mixing. 

To  visualize  the  capture  zone  and  the  shadow  zone  a  head  difference  of  4.2mm  was  applied  to 
the  0.50mm  aperture  top  fracture  resulting  in  a  O.lOml/s  volumetric  flow  rate  through  the  241mm 
wide  fracture  and  a  velocity  of  0.83mm/s.  Dye  was  distributed  uniformly  within  the  constant 
head  reservoir  and  tracked  as  it  traversed  the  fracture  and  entered  the  borehole.  Initial  runs  did 
not  control  for  temperature  or  density  differences  between  the  fracture  and  the  borehole.  These 
temperature  and  density  contrasts  allow  the  characteristics  of  density  mixing  to  be  observed. 


4,3  Analytical  and  Numerical  Modeling 

4. 3. 1  Analytical  Modeling  Approach  and  Methods 

Analytical  models  of  wellbore  hydrodynamics  provide  a  rapid  and  flexible  method  of 
investigating  fundamental  aspects  of  wellbore  flow  during  pumping.  Analytical  models  also 
allow  the  results  to  be  generalized  in  order  to  provide  design  criteria  for  pumping  and  sampling 
(Martin-Hayden  2000b).  Previous  models  assumed  inviscid  flow,  i.e.,  Darcian  flow,  within  the 
wellbore.  These  models  provide  estimates  of  travel  time  within  the  wellbore  and  relationships 
between  the  groundwater  concentration  distributions  and  concentration  variation  during 
pumping.  In  this  study  more  complex  analytical  and  numerical  models  are  developed  to  account 
for  Poisueille  (viscous  pipe)  flow  and  investigate  the  influences  of  viscosity  on  flow  paths,  travel 
time  and  concentration  distributions.  The  velocity  distributions  within  the  screen  section  of  the 
well  are  then  compared  to  various  mechanisms  of  wellbore  mixing,  i.e.,  thermal  and  solutal 
convection  as  well  as  advection-diffusion/dispersion. 

The  analytical  modeling  began  with  investigations  of  wellbore  flow  during  pumping.  The 
distribution  of  groundwater  entering  the  screen  was  modeled  as  an  advective  front  forming  an 
annulus  of  groundwater  displacing  a  shrinking  cylinder  of  pre-pumping  well  water.  If  the 
groundwater  and  well  water  are  heterogeneous  in  concentration,  and  the  well  is  unmixed  or 
partially  mixed  during  pumping,  the  distribution  of  groundwater  and  well  water  is  used  to 
characterize  the  heterogeneity  of  well  water  concentration  during  pumping.  The  travel-time 
calculations  show  how  the  concentrations  at  the  pump  intake  will  change  with  time  (i.e., 
“pumping  curves”)  and  can  be  used  to  characterize  variation  of  the  pumped  concentration. 
These  pumping  curves  were  generated  for  uniform  groundwater  concentrations  and 
heterogeneous  concentrations  with  the  pump  intake  at  various  positions  within  the  screened 
interval  in  order  to  simulate  a  variety  of  conditions  during  pumping. 

When  a  screened  monitoring  well  is  not  pumped  (between  sampling  episodes,  if  sampled  with  a 
passive  method  or  monitored  with  an  in  situ  sensor)  it  acts  as  a  highly  conductive  conduit  that  is 
very  sensitive  to  wellbore  mixing  and  vertical  flow.  The  concentrations  carried  to  the  well 
depend  on  the  surrounding  flow  regime  established  by  the  interaction  of  the  well  with  the 
ambient  hydraulic  gradients,  vertical  as  well  as  horizontal  gradients.  The  redistribution  of 
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concentration  within  the  well  depends  on  the  vertieal  flow  in  addition  to  the  meehanisms  driving 
mixing.  Vertieal  flow  and  various  factors  driving  convection  will  work  either  in  coneert  or  in 
opposition  to  redistribute  groundwater  within  the  well.  In  this  study,  modeling  of  the  eonduetive 
well  bore  and  the  surrounding  formation  showed  the  flow  regime  carrying  groundwater  to  the 
well,  vertical  flow  within  the  well  and  the  distribution  of  well  water  discharge  creating  the 
shadow  of  altered  and  redistributed  water. 

4.3.2  Numerical  Modeling  Approach 

Numerieal  modeling  studies  have  been  condueted  to  examine  how  solute  distribution  in  a  well  is 
influenced  by  low-flow  pumping  and  by  thermal  and  solutal  eonvection  effeets  and  to  examine 
the  possible  influenee  of  well-sereen  geometry  on  in-well  solute  transport.  Many  previous 
studies  of  groundwater  flow  and  solute  transport  that  have  included  in-well  transport  effeets  have 
investigated  flow  under  pumping  conditions  in  which  flow  is  dominated  by  the  pumping-indueed 
foreed  eonveetion  (Reilly  and  Gibs,  1993;  Chiang  et  ah,  1995;  Elei  et  al.  2001;  Varljen,  2006). 
Under  those  conditions,  flow  in  the  well  can  be  approximated  by  solutions  that  use  Darey’s  law 
to  deseribe  conservation  of  momentum.  That  assumption  may  be  invalid,  however,  when  in-well 
veloeities  due  to  pumping  are  of  the  same  order  of  magnitude  as  flow  rates  imposed  by  other 
forees,  sueh  as  thermal  or  solutal  convection.  Development  and  use  of  in  situ  monitoring  deviees 
would  normally  require  monitoring  in-well  eoneentrations  without  purging  or  pumping.  Under 
those  conditions  thermal  or  solutal  eonveetion  effeets  may  exert  significant  influence  on  in-well 
solute  transport.  As  such,  an  accurate  understanding  of  that  potential  influence  is  required  to 
understand  how  in-well  solute  monitoring  would  reflect  solute  distribution  in  the  surrounding 
formation.  Navier-Stokes  equations  for  fluid  flow  are  needed  to  deseribe  flow  and  transport  in 
the  well  under  these  eonditions.  However,  even  where  flow  is  dominated  by  pumping,  there  are 
eireumstanees  where  application  of  Darcy’s  law  may  yield  inaeeurate  results.  One  sueh  example 
is  the  displacement  of  initial  well  water  during  well  purging,  for  which  Darcy’s  law  provides  an 
ineorrect  assumed  stream-tube  flow  distribution  and,  therefore,  an  inaeeurate  deseription  of 
purging  behavior  during  pumping. 

To  more  accurately  deseribe  in-well  flow  and  transport  behavior  and  to  incorporate  potential 
buoyaney  effeets,  we  simulated  flow  and  transport  using  the  3D  finite  element  analysis  program 
COMSOL  Multiphysies.  The  Multiphysies  paekage  allows  solution  of  fully  eoupled  fluid  flow, 
solute  transport  and  thermal  energy  transport  on  a  finite-element  domain,  with  fluid  flow 
governed  by  Darcy’s  law,  the  Brinkman  equations  or  Navier-Stokes  equations  for  weakly 
eompressible  flows.  The  software  provides  convenient  pre-  and  post-processing  tools,  and  for 
problems  involving  solution  of  different  governing  equations,  the  solutions  are  eoupled  at  the 
boundaries  between  the  combined  solution  domains. 

4.3.2  Numerical  Modeling  Methods 

It  is  eommonly  assumed  that  laminar  flow  through  a  well  forees  the  solute  distribution  inside  a 
well  to  mirror  that  outside  the  well.  This  would  suggest  that  in  cases  where  solute  is  confined  to 
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certain  strata  in  a  hydrogeologic  profile,  that  same  vertical  distribution  would  be  refieeted  in  a 
monitoring  well.  To  test  that  conceptual  model  of  flow  through  a  well,  Britt  (2005)  eonstructed  a 
2D  sand  tank  that  ineorporated  a  vertical  well  section  in  a  horizontal  flow-dominated  system 
(Figure  6).  The  well  consisted  of  a  lO-em  diameter  included  a  eoarse  sand  paek  and  divided  two 
porous  medium  domains  in  whieh  horizontal  flow  was  controlled  by  eonstant  head  tanks  on 
either  end  of  the  apparatus.  A  dye,  mixed  with  water  to  yield  a  solution  density  -l-/-0.00001x  that 
of  water,  was  injected  just  upstream  of  the  well  to  examine  solution  behavior  as  the  dye  migrated 
laterally  toward  and  into  the  well.  In  all  experiments  the  dye  was  observed  to  migrate  vertieally 
in  the  well  much  faster  than  horizontally.  While  the  horizontal  rate  of  dye  movement  was  well 
approximated  by  the  volumetric  flux,  the  vertieal  rate  was  on  the  order  of  30  times  faster.  INL 
condueted  experiments  that  essentially  replicated  Britt’s  apparatus,  at  similar  horizontal 
velocities  and  up  to  an  order  of  magnitude  greater  than  the  2005  study. 


P3070017@  3/07/03  817a 
test  start  420p  on  3/06/03 
calc,  dye  density:  0.999986  g/cc  ] 
7:00  hrs.  since  dye  emerged 
I  seepage  velocity:  0.47  ft/day 


Figure  6.  The  quasi-2D  horizontal  flow  to  a  well  apparatus  of  Britt  (2005),  used  to  examine  solute  transport  through 
a  well  under  ambient  horizontal  flow.  Far  left  and  right  columns  provided  constant  head  boundaries  while  sealed 
upper  and  lower  surfaces  result  in  a  conflned  hydraulic  system.  Dye  was  injected  at  the  marked  port  and  migrated  to 
the  well  at  a  velocity  of -0.47  ft  per  day,  with  limited  dispersion  approximately  denoted  by  indicators  on  the  edge  of 
the  visible  plume  in  the  sand.  Dye  distribution  illustrates  rapid  vertical  movement 

Britt’s  (2005)  experiments  demonstrated  that  slight  density  differenees  did  yield  different 
vertieal  distribution  of  solutes  in  the  well.  To  estimate  the  sensitivity  of  a  well  to  density-driven 
mixing  effeets,  we  use  dimensionless  parameters  that  compare  the  forees  driving  buoyancy- 
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driven  mixing  to  those  forees  that  attenuate  that  effect.  Effects  that  attenuate  flow  due  to  a  net 
buoyant  force  include  viscosity  and  diffusivity,  and  the  Rayleigh  number  is  the  dimensionless 
parameter  that  expresses  the  balance  between  those  effects.  The  form  of  the  Rayleigh  number 
and  the  value  at  which  natural  convection  begins  varies  with  the  geometry  of  a  particular  system. 
For  flow  in  a  vertical  tube  separating  two  solutions  of  differing  density  (Figure  ),  the  Rayleigh 
number,  Ra,  is 


where  g  is  gravitational 
acceleration,  R  is  radius  of  the 
tube,  p,  is  kinematic  viscosity, 

D  is  the  diffusion  coefficient,  P 
is  the  fractional  volumetric 
expansion  coefficient,  and 
dc/dz  is  the  concentration 
maintained  by  the  reservoirs  on 
each  end  (Cussler,1997). 

Taylor  (1954)  determined  that 
the  critical  Rayleigh  number 
for  the  onset  of  natural 
convection  for  the 

configuration  shown  in  Figure 
is  68.  Thus,  for  a  10-cm 
diameter  tube,  the  critical  concentration  gradient  for  solute  driven  convection  is  approximately 
0.4  pg/(F  m).  In  dilute  solutions,  the  volumetric  expansion  coefficient  associated  with  slight 
increases  in  solute  concentration  is  approximately  lE-3  F/kg,  so  the  relative  density  difference 
represented  by  this  concentration  gradient  is  approximately  4E-10  m'^ 

Thermal  gradients  can,  of  course,  also  drive  convective  flow,  but  because  heat  can  also  diffuse 
into  the  well  and  surrounding  medium  the  critical  Rayleigh  number  depends  on  the  ratio  of  the 
thermal  conductivity  of  the  fluid  in  a  tube  to  that  of  the  surrounding  medium.  Zhukhovitskii  and 
Gershuni  (1976)  show  that  the  critical  Rayleigh  number  for  thermal  convection  in  a  tube  can  be 
calculated  from  the  expression: 

Ra  =  +  103;.)  -  ^3(2567  +  14794;  +  26927A=)J, 

where  X  is  the  ratio  of  the  thermal  conductivities  inside  and  outside  the  tube.  When  the  thermal 
conductivity  of  the  fluid  is  much  greater  than  the  surrounding  medium,  this  equation  yields  a 
critical  Rayleigh  number  of  68,  identical  to  the  solutal  convection  case.  For  conditions  more 
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typical  of  the  saturated  zone,  the  critical  Rayleigh  number  is  approximately  149  (Berthold  and 
Bomer,  2008).  Using  the  latter  number  to  estimate  when  natural  convection  would  dominate 
transport  by  diffusion  in,  again,  a  5-cm  diameter  well,  we  see  that  a  thermal  gradient  of  0.0035 
K/m  would  be  sufficient,  which  is  approximately  an  order  of  magnitude  lower  than  typical 
geothermal  gradients.  The  fractional  expansion  coefficient  for  water  at  10°C  is  approximately 
8.8E-5  K'\  so  the  critical  gradient  reflects  a  relative  density  difference  of  approximately  3E-7  m' 
\  This  is  approximately  three  orders  of  magnitude  larger  than  necessary  to  produce  solutal 
convection  because  the  dissipative  effect  of  thermal  diffusivity  is  approximately  that  much  larger 
than  typical  solutal  diffusion  coefficients  (Eigure  8). 

Because  the  Rayleigh  number  is  proportional  to  the  fourth  power  of  tube  radius,  the  critical 
gradient  describing  the  onset  of  convection  has  the  non-linear  dependence  on  radius  shown  in 
Figure  8. 
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Figure  8.  Critical  gradient  for  the  onset  of  convection  in  a  tube,  for  thermal  and  solutal  convection,  as  a  function  of 
tube  radius. 

The  Rayleigh  number  analysis  defines  when  free  convection  will  dominate  diffusion  in  thermal 
energy  or  solute  transport  but  it  does  not  predict  the  scale  of  the  convective  cells  that  develop  or 
the  rate  of  mixing  induced  by  convection  (Gretener,  1967).  This  information  is  needed  to  predict 
the  effect  of  natural  convection  on  solute  distribution  in  a  well,  particularly  where  forced 
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convection  in  the  horizontal  or  vertical  direction  may  also  influence  solute  movement. 
Parametric  analysis  of  geometrically  described  flows  adds  some  insight  for  particular  scenarios. 
For  flow  in  a  cylinder,  for  example,  the  flrst  convective  flow  regime  expected  to  occur  is  flow 
down  one  half  of  the  cylinder  with  a  compensating  upward  flow  along  the  opposite  half  of  the 
cylinder  (Donaldson,  1961).  Axisymetric  flow,  with  one  flow  direction  along  the  center  of  the 
cylinder  and  the  oppositely  directed  flow  in  an  annular  ring  surrounding  it,  does  not  occur  until 
higher  Rayleigh  numbers.  To  examine  the  details  of  convective  flow  behavior,  we  use  numerical 
simulation,  solving  the  equations  for  open  water  flow  (Navier-Stokes)  in  domains  representing 
wells  and  applying  Darcy’s  law  in  porous  media.  Heat  transport  and  solute  transport  are 
calculated  by  coupling  the  advection-dispersion  equation  to  the  local  flow  equation,  which  may 
vary  between  domains. 

Numerical  simulations  of  solute  transport  were  conducted  using  several  different  domains, 
depending  on  the  problem  of  interest.  To  examine  solute  transport  and  heat  transport  in  an 
experimental  apparatus  like  that  of  Britt  (2005),  we  developed  models  that  essentially 
reproduced  the  2-D  geometry  and  boundary  conditions  of  that  experiment  (Figures  9,  lOA),  or 
used  a  similar  geometry  with 
slightly  taller  aspect  (e.g. 

Figure  lOB),  to  examine 
vertical  movements  of  greater 
extent.  In  these  and  other 
models,  flow  in  the  well 
section  was  described  using 
Navier-Stokes  equations,  while 
flow  in  the  surrounding  porous 
media  domains  was  described 
using  Darcy’s  law.  Solute 
transport  and  heat  transport 
were  described  using 
advection-dispersion  models 
and  -  in  cases  without  forced 
convection  -  density-based 
buoyant  forces  included  using  a 
Boussinesq  approximation. 

To  examine  solute  transport  under  pumping  conditions,  we  simulated  flow  and  transport  in  a  1-m 
by  10-cm  radius  well,  surrounded  by  a  10-cm  zone  of  porous  media  (Figure  11  A),  using  an 
axisymmetric  domain.  Pumping  in  that  model  was  represented  as  a  specifled  flux,  concentration 
outflow  boundary  along  the  bottom  of  the  well  section. 
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Figure  9.  Schematic  of  the  finite  element  model  used  to  simulate  Britt’s 
(2005)  experiment  that  displayed  strong  vertical  mixing  with  relative 
density  gradients  less  than  lE-5. 
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The  effect  of  time-varying  temperature  gradients  on  flow  in  a  well  were  also  investigated  with  a 
1-m  tall  well  section  in  an  axi-symmetric  model  domain  (Figure  IIB),  but  with  an  extended 
aquifer  surrounding  the  well  to  allow  more  accurate  description  of  the  transient  heat  exchange 
between  the  well  and  the  aquifer. 


Figure  10.  Model  domain  geometry  (distances  in  meters)  and  boundary  conditions  for  simulations  of  Britt’s  2005 
laboratory  experiment  (A)  and  similar  scenarios  with  greater  vertical  extent  (B). 
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Figure  1 1 .  Model  domain  geometry  (distances  in  meters)  and  boundary  conditions  for  simulations  of  solute  transport 
in  a  well  with  extraction  from  any  point  along  the  screen  (A)  and  free  convection  behavior  in  a  well  under  constant 
or  time-varying  temperature  gradients  induced  by  seasonal  temperature  fluctuations  (B). 


4,4  Field  Experiments 

Field  experimentation  is  a  key  factor  in  establishing  the  predictive  capability  of  the  physical  and 
numerical  modeling.  Many  aspects  of  the  experimental  laboratory  work  require  field  validation. 
Additional  field  experimentation  to  test  in-well  mixing  effects,  as  well  as  effects  of  purge  rate, 
purge  rate  change,  density  effects  and  other  parameters  were  needed  to  establish  the  dynamics  of 
variability  associated  simply  with  the  way  traditional  samples  are  collected. 

A  field  site  for  the  study  was  selected  and  approval  for  its  use  was  obtained  from  the  private  site 
owner  and  the  local  oversight  agency  (the  California  Department  of  Toxic  Substances  Control). 
The  field  site  is  the  former  Angeles  Chemical  Facility  located  at  8915  Sorenson  Avenue  in  Santa 
Fe  Springs,  California  (Figure  12).  The  site  has  a  wide  range  of  volatile  organic  compound 
contamination,  including  petroleum  hydrocarbons,  chlorinated  solvents,  ketones,  and  1,4- 
dioxane.  Contaminant  concentrations  range  from  below  detection  level  to  NAPL-indicating 
concentrations.  Hydrogeologicahy,  the  site  has  a  range  of  flow  regimes  from  nearly  stagnant, 
poor-yield,  and  highly  heterogeneous,  to  high  flow,  high-yield,  and  low-heterogeneity.  Wells 
include  short  to  moderate  length  screens  5-15  feet,  of  2-inch  and  4-inch  construction.  The  site  is 
less  than  2  acres  in  total  area.  A  Remedial  Investigation  Report  prepared  for  the  site  by  the 
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California  Department  of  Toxic  Substances  Control  was  available  to  guide  the  site  work.  Details 
of  the  site  geology  and  investigation  derived  data  are  included  in  that  report. 


Historically,  several  individual  wells  at  the  site  have  been  sampled  at  intervals  to  identify 
contaminant  stratification.  Many  were  suitable  candidates  for  active  and  passive  sampling 
approaches  to  illustrate  laboratory  experimental  findings.  During  the  first  field  season  we 
conducted  the  following  tests: 

•  Multilevel  passive  sampling.  Multilevel  sampling  utilizing  existing  site  wells  was 
conducted  using  two  approaches:  1)  open  well  sampling  with  several  stacked  passive 
sampling  devices  deployed  in  each  open  well  to  assess  ambient  stratification;  2)  isolated 
zone  sampling  utilizing  stacked  passive  samplers  deployed  in  each  well  with  baffie 
isolation  devices  to  limit  ambient  mixing  between  zones. 

•  Tracer  Tests.  After  the  passive  samplers  were  removed  from  the  existing  test  wells, 
several  wells  were  outfitted  with  sensor  devices  for  tracer  testing.  The  tracer  (deionized 
water)  was  placed  in  the  wells  to  measure  flow  through  via  tracer  dissipation. 

•  Well  Installation.  Based  on  passive  sampling  results  and  tracer  tests,  new  wells  were 
installed  in  locations  where  contaminant  stratification  was  identified.  Short  screen  wells 
along  with  tube  well  ports  within  the  well  filter  pack  were  incorporated  into  the  well 
installation.  Sampling  to  explore  the  relationship  between  inside  and  outside  the  well 
during  ambient  flow  and  during  pumping  was  used  to  illustrate  small-scale  vertical 
heterogeneity. 

•  Purge  Dynamics  Testing.  Low  flow  purging  and  volume-based  purging  comparisons 
included  protocols  to  test  how  pumping  rate,  consistency  of  flow  rate,  changes  in  flow 
rate  prior  to  sampling,  and  duration  of  pumping  affects  pump  discharge  concentrations, 
relative  to  concentrations  measured  in  the  well.  Physical  testing  of  flow  dynamics  within 
individual  test  wells  was  accomplished  through  high  resolution  temperature,  ORP,  and 
conductivity  tools,  and  other  physical  methods  to  define  physical  conditions  and 
dynamics  downhole.  These  physical  experiments  were  done  in  conjunction  with 
chemical  testing  to  illustrate  flow  patterns. 
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•  Bottle  Filling  Test.  In  Addition  to  the  passive  sampling  and  purge  curve  sampling 
analysis,  a  single  outdoor  field  test  of  methods  that  are  commonly  used  to  fill  VO  A  vials 
was  conducted.  This  test  was  conducted  to  assess  the  degree  to  which  field  sample 
handling  affects  analytical  results. 


4.4.1  Multilevel  Passive  Sampling 

Eight  existing  wells  at  the  site  were  sampled  at  multiple  levels  within  the  saturated  screen 
interval.  The  intent  of  the  initial  sampling  episode  was  to  determine  current  levels  of 
contaminant  stratification  in  local  geographic  areas  of  the  selected  wells.  Stratification 
information  was  used  to  select  locations  for  new  well  installations  required  to  accommodate 
equipment  for  the  field  program.  Two  modes  of  passive  sampling  were  employed  for  the  initial 
sampling  phase.  One  mode  included  passive  sampling  in  the  open  well,  while  the  second  mode 
included  zone  isolation  between  sampling  devices.  Snap  Samplers  were  selected  for  the  passive 
sampling.  Of  the  available  passive  samplers,  the  Snap  Sampler  is  the  only  device  that  collects  an 
undisturbed  whole-water  sample  of  the  ambient  flow-through  present  in  the  well  at  the  time  of 
sampling.  Multiple-interval  samples  are  achievable,  and  installation  of  zone  separation  baffles  is 
relatively  easy  to  accomplish  using  this  device  (Figure  13).  The  following  wells  were  selected 
for  preliminary  interval  monitoring:  MW-9,  MW-14,  MW-15,  MW-16,  MW-17,  MW-20,  MW- 
21,  and  MW-24.  Figure  14  shows  locations  of  these  wells. 


Figure  13.  In  well  flow  inhibitor  used  by 
ProHydro,  Inc.  to  limit  vertical  solute  migration 
in  a  well. 
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Figure  14.  Well  Location  Map 

A  total  of  26  sets  of  samplers  were  deployed  in  the  selected  wells.  Each  well  was  outfitted  with 
2  to  4  sets  of  Snap  Samplers  activated  with  pneumatic  trigger  systems.  Pneumatic  trigger  lines 
were  fixed  at  each  well  head  to  assure  they  remain  in  a  fixed  position  during  sampler 
equilibration.  Trigger  lines  were  measured/length-corrected  in  the  field  to  accommodate 
discrepancies  from  available  well  construction  information.  Deployment  was  completed 
according  the  depth  intervals  shown  in  Table  1.  Because  up  to  four  zones  were  monitored  in  10- 
foot  well  screens,  and  two  zones  monitored  in  the  5  foot  saturated  screens,  a  maximum  of  only 
two  Snap  Samplers  were  deployed  with  the  actuator  in  any  single  2.5-ft  interval.  Most  samples 
consisted  of  (2)  40  ml  VGA  vial  Snap  Samplers.  For  duplicate  samples.  Matrix  Spikes,  and 
Matrix  Spike  Duplicate  samples,  samples  were  collected  using  either  125ml  or  350ml  Snap 
Sampler  bottles.  For  these  QA/QC  samples,  each  sample  was  poured  into  standard  laboratory- 
supplied  VGA  vials.  By  collecting  the  QA/QC  samples  in  this  manner,  the  replicates  are  true 
replicates  rather  than  a  combination  of  Snap  VGAs  and  lab  VGAs.  This  also  eliminated 
concerns  with  spatial  variability  or  collection  and  handling  methods.  In  this  initial  installation. 
Snap  Samplers  were  deployed  with  isolation  baffles  separating  each  sampling  zone. 
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Table  1,  Snap  Sampler  Deployment  Positions 


Well 

MW-9 

MW-14 

MW-15 

MW-16 

MW-17 

MW-20 

MW-21 

MW-24 

TOC  Elevation 

149.16 

149.71 

150.6 

148.32 

149.03 

149.14 

150.02 

149.9 

screen  interval 
(bgs) 

30.5-45.5 

55-65 

54-64 

29-46 

56-66 

57-67 

53-63 

67-77 

Well  diameter 
(in) 

4 

2 

2 

2 

2 

2 

2 

4 

Depth  to  Water 
(ft)  03/09 

40.31 

54.2 

55.64 

39.9 

52.2 

52.95 

54.16 

53.94 

Sample  zone 

Zone  1 

41 

56 

57 

41 

57 

58 

55.5 

68 

Zone  2 

43.5 

58.5 

59.5 

43.5 

59.5 

60.5 

58 

70.5 

Zone  3 

61 

62 

62 

63 

61.5 

73 

Zone  4 

63.5 

64.5 

65.5 

75.5 

Following  an  approximate  two-week  equilibration  period,  Snap  Samplers  were  triggered  and 
retrieved  (1st  week  of  June  2009).  Upon  eompletion  of  initial  sampling  aetivities  at  eaeh  well, 
each  isolation  baffle  was  removed  prior  to  redeployment.  The  purpose  of  the  second  sampling 
event,  and  for  removing  the  baffles,  is  to  test  the  difference  in  ambient  stratification  the  open 
well.  Samplers  were  deployed  and  again  allowed  to  equilibrate  within  the  wells  for 
approximately  two  weeks.  After  the  second  sampling  event  equilibration  period.  Snap  Samplers 
were  triggered,  retrieved  and  prepared  for  laboratory  submittal. 

Later  in  the  field  program  (2011),  a  second  multilevel  field  testing  program  was  conducted.  The 
second  phase  included  deployment  of  polyethylene  diffusion  bag  samplers  at  1  foot  intervals  in 
two  wells.  The  two  wells  were  installed  with  multilevel  sampling  ports  (“tube  wells”)  described 
in  more  detail  below.  The  PDBs  were  isolated  from  each  other  using  plastic  baffle  devices 
designed  to  limit  in  well  mixing  effects.  The  purpose  of  these  test  was  to  compare  the  outside 
tube  well  samples  to  isolated  passive  devices  located  inside  the  well. 

4.4.2  Tracer  Tests 

After  removal  of  all  of  the  sampling  equipment  from  the  first  passive  sampler  deployment  in 
2009,  the  eight  existing  test  wells  were  tested  for  natural  flushing  using  a  tracer.  Downhole 
logging  sensors  were  deployed  in  two  depths  (in  each  well)  to  collect  tracer  dissipation  data. 
Table  2  provides  the  sensor  deployment  depths  for  all  the  sensors.  The  tracer  deionized  water 
and  was  emplaced  using  a  tremie  to  introduce  minimal  agitation,  yet  promote  thorough  mixing  of 
the  tracer  in  the  well  screen  column.  Minimally  sufficient  tracer  was  applied  to  allow 
measurement  without  significant  disruption  of  either  the  flow  regime  or  long  term  aquifer 
chemistry.  Recovery  of  well  chemistry  to  the  pre -tracer  condition  was  recorded  with  the  logging 
sensors.  Sensors  measured  temperature,  pH,  oxidation-reduction  potential,  electrical 
conductivity  and  chloride.  Data  was  collected  for  up  to  two  weeks,  depending  on  field  schedule. 
Measurements  were  logged  at  15  minute  intervals  over  the  course  of  deployment  period  to  allow 
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sufficiently  frequent  measurement  of  tracer  dissipation.  Flushing  rate  was  determined  from  this 
tracer  approach. 


Table  2,  Sensor  Deployment  Positions 


Well 

MW-9 

MW-14 

MW-15 

MW-16 

MW-17 

MW-20 

MW-21 

MW-24 

TOC  Elevation 

149.16 

149.71 

150.6 

148.32 

149.03 

149.14 

150.02 

149.9 

screen  interval 
(bgs) 

30.5- 

45.5 

55-65 

54-64 

29-46 

56-66 

57-67 

53-63 

67-77 

Well  diameter  (in) 

4 

2 

2 

2 

2 

2 

2 

4 

Depth  to  Water  (ft) 
03/09 

40.31 

54.2 

55.64 

39.9 

52.2 

52.95 

54.16 

53.94 

Sensor  zone 

Zone  1 

41 

57.5 

57 

41 

58.5 

60.5 

56 

70.5 

Zone  2 

43.5 

62.5 

61 

43.5 

63.5 

65.5 

60 

75.5 

4.4.3  Installation  of  New  Wells 

The  drilling  program  was  designed  to  accommodate  sampling  requirements  of  the  active  pump¬ 
sampling  phase  of  the  investigation.  Most  of  the  on-site  wells  are  2-inch  diameter  and  would  not 
accommodate  all  of  the  equipment  required  to  run  the  purging  dynamics  tests.  Additionally,  the 
testing  program  required  sampling  ports  located  outside  the  wells,  in  the  annular  space  between 
the  well  and  the  borehole  wall.  Therefore  three  new  wells  were  installed  to  satisfy  those 
requirements. 

The  new  wells  are  4-inches  in  diameter,  with  multilevel  sampling  ports  installed  outside  the  well 
casings  (Table  3).  New  well  locations  and  specific  design  criteria  were  based  on  results  of  the 
initial  well  sampling  and  tracer  testing  described  above.  Figure  15  illustrates  the  construction 
detail.  Wells  were  installed  by  standard  methods,  and  developed  thoroughly  by  surging,  bailing 
and  pumping  to  achieve  high  quality  sample  access  points.  During  development,  sustainable 
flow  rates  for  the  wells  were  determined. 

Table  3,  New  Well  Construction  Details 


Well 

MW-27 

MW-28 

MW-29 

Depth  (ft) 

65 

62.5 

46 

Well  Diameter  (in) 

4 

4 

4 

Diameter  of  hole  (in) 

18/12 

18/12 

12 

Screen  interval  (ft) 

55-65 

52.5  -62.5 

41-46 

Tnhe  well  depths  (ft) 

55.5- 56.5;  58-59; 

60.5- 61.5;  63-64 

53-54;  55.5-56.5; 
58-59;  60.5-61.5 

42,  45 

Material  screened 

medium-coarse 
clean  sand 

medium-coarse 
clean  sand 

silty  sand 
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Flush-threaded  casing  Joints 


10  ft  4"  sch40  PVC  0.020  slot 


RMC2/12  Sand,  Ift  above  TOS 


Figure  15.  Construction  Detail  for  New  Wells 

After  the  wells  were  installed  and  developed,  the  wells  were  instrumented  with  sensors  and 
outfitted  with  Grundfos  RediFlo2  pump  or  bladder  pump  equipment  and  Snap  Samplers  in 
antieipation  of  purge  sampling  work.  The  multi-parameter  sensors  were  eapable  of  measuring 
and  reeording  pressure  (depth  of  submergence),  electrical  conductivity,  temperature,  oxidation- 
reduction  potential,  chloride  ion  and  pH.  The  RediFlo2  is  an  electric  submersible  pump  capable 
of  adjustable  flow  rates  as  low  as  approximately  lOOml/minute  up  to  several  liters  per  minute. 
The  bladder  pump  used  in  MW-29  (a  low  recharge  well)  was  capable  of  flow  rates  as  low  as 
50ml  per  minute  up  to  approximately  250ml  per  minute.  Equipment  installed  in  the  wells  was 
left  to  stabilize  for  one  week  to  one  month  in  advance  of  sampling.  The  stabilization  period  is 
designed  to  allow  the  well  itself,  and  the  equipment  installed  to  chemically  equilibrate  with  the 
surrounding  formation  water.  Previous  research  (Parker  and  Ranney  1997,  1998,  2000,  2003) 
shows  that  new  plastics  can  sorb  contaminants  and  thus  could  bias  results  if  not  equilibrated. 
Pre-deployment  of  the  sensors,  pumps,  samplers  will  limit  potential  bias  from  new  materials  in 
the  equipment.  Figure  16  graphically  represents  the  arrangement  of  equipment  to  be  deployed  in 
the  test  wells  for  the  purge  dynamics  testing. 
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Figure  16.  Downhole  Equipment  Illustration  and  Photo. 


The  tube  sampling  ports  were  placed  on  the  upgradient  side  of  the  well  during  installation.  By 
installing  the  wells  and  equipment  well  in  advance  of  sampling,  natural  flow-through  is  allowed 
to  establish  a  dynamic  equilibrium  between  the  equipment  and  flowing  formation  water. 

4.4.4  Purge  Dynamics  Testing 

Several  tests  were  completed  during  the  purge  testing  phase.  Low  flow  purging  (defined  as  any 
flow  rate  less  than  approximately  1000ml  per  minutes  that  generates  less  than  0.3ft  of 
drawdown)  was  conducted,  with  multiple  VOC  samples  collected  along  the  “purge”  curve. 
These  tests  were  designed  to  test  how  flow  rate,  flow  rate  changes,  and  variable  contaminant 
concentrations  along  the  well  screen  influencing  aggregate  concentrations  pumped  from  the  well, 
and  other  changes  observed  during  the  course  of  purging.  Multiple  in-well  sensors  measured 
water-quality  parameters  (pH,  ORP,  EC,  Cl)  of  the  inflowing  water,  chemistry  as  well  as  these 
parameters  were  measured  in  the  pumped  discharge  water  as  well.  Before  purging  and  toward 
the  end  of  a  purge  cycle,  two  sets  of  Snap  Samplers  were  triggered  to  compare  water  flowing  to 
the  pump  from  above  and  below  the  pump  intake,  before  and  after  purging.  Tube  wells  located 
at  depths  along  the  screen  interval,  just  outside  the  well  screen  were  used  to  collect  samples  of 
water  moving  to  the  well  at  the  end  of  the  purge  cycle.  Figure  17  illustrates  flow  paths  of  water 
entering  the  well  during  pumping. 

Varljen,  et  al.  (2006)  simulated  steady  state  pumping  scenarios  numerically  at  steady  state,  with 
similar  graphic  representations.  Varljen  posited  that  since  steady-state  inflow  came  from  the 
whole  screen  zone,  the  discharge  is  basically  by  definition  a  flow  weighted  average.  Martin- 
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Hayden  (2000a)  provided  an  analytie  solution  for  inflowing  eontaminant  coneentrations  over 
time,  whieh  made  a  case  for  strong  influence  on  pump  discharge  by  inflowing  contaminant 
position  relative  to  the  pump  position  and  to  time  of  purging.  The  experiments  in  this  project, 
with  instrumented  sampling  wells,  were  designed  to  help  clarify  how  the  variability  of  various 
factors  (aquifer  permeability,  contaminant  concentrations,  pump  location,  extraction  rate,  and 
time)  influence  what  pump  discharge  “represents.” 


Illustration  only:  Flow  path  stait  positions 
are  NOT equivalentto  traveltiine 


Figure  17.  Flow  Paths  During  Active  Pumping 


As  described  above,  the  3  new  test  wells  were  instrumented  with  sensors.  Snap  Samplers  and  a 
purge  pump  several  days  to  a  month  in  advance  of  purge  sampling.  The  Snap  Samplers  were 
deployed  just  above  and  below  the  pump  intake  position.  Sensors  were  deployed  at  the  same  2 
positions  as  the  Snap  Samplers  and  at  third  interval  (above)  and  fourth  interval  (below)  the  Snap 
Samplers.  The  sensor  positions  were  approximately  the  same  as  the  tube  well  sampling  ports. 
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This  positioning  was  designed  to  collect  data  on  water  flow/chemistry  ahead  of  purge  sampling 
activities,  and  to  document  changes  as  purge  sampling  progresses. 

Four-inch  well(s)  have  a  nominal  inside  diameter  of  10cm  and  a  length  of  305cm.  One  fully 
saturated  well  casing  volume  is  therefore: 

IWV  =  (lOcm/2)^  X  3.14  X  305cm 

IWV  =  23.9  X  10^  cm^,  or  23.9  liters 

For  2-inch  (5cm)  wells: 

IWV  =  (5cm/2)^  X  3.14  X  305cm 

IWV  =  6.25  X  10^  cm^,  or  6.0  liters 

Purge  sample  protocol 

•  Start  sensor  parameter  logging  at  1  minute  intervals  approximately  1  hour  before  pump 
start  (continue  for  duration  of  all  testing  activities). 

•  Activate  first  set  of  Snap  Samplers  to  close. 

•  Start  pumping  lOO-lOOOmL/minute  as  determined  during  well  development  for  steady 
drawdown. 

•  Collect  VOC  samples  from  the  pump  discharge  at  a  rate  of  approximately  6  samples  per 
well  volume,  or  one  per  4  liters  for  5  well  volumes. 

•  For  wells  with  a  flow  rate  change  test,  increase  purge  rate  to  a  rate  that  will  increase 
drawdown  approximately  2-3x  that  of  the  lower  rate,  as  determined  during  well 
development. 

•  For  wells  with  a  flow  rate  change  test,  reduce  flow  back  to  initial  rate,  to  allow  water 
level  to  rise  during  purging;  collect  VOC  samples  every  1  minute  for  10  minutes, 
followed  by  samples  every  5  minutes  for  20  more  minutes. 

•  Collect  annular  space  sampling  port  samples  (“tube  well  samples”)  during  the  increased 
flow  rate  purging  episode. 

•  Snap  the  second  set  of  Snap  Samples,  turn  off  pump,  &  recover  Snap  Samplers  from  the 
well. 

Table  4  includes  a  listing  of  the  wells,  flow  rates,  times  of  purging,  volumes  purged  and  number 
and  type  of  samples  collected  during  the  tests  conducted  during  the  trial.  Appendix  A  contains 
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the  details  of  the  sampling  events  and  samples  collected  during  the  purge  intervals.  Analytical 
results  will  be  described  in  greater  detail  in  the  Results  section. 

Table  4,  Samples  Collected  During  Purge  Testing 


Well 

(diameter) 

Date 

Flow 

Rate 

mFmin 

Minutes 

purged 

Total 

Volume/ 

well 

volumes 

purged 

Number 

of 

Samples 

Snap 

Samples 

Tube 

Samples 

Duplicate 

and/or 

MS/MSD 

MW-24  (4”) 

8/24/09 

840 

158 

213/9 

42 

4 

N/A 

n 

MW-28  (4”) 

8/25/09 

900 

98 

92/4 

33 

4 

2 

y(2) 

MW-29  (4”) 

9/1/09 

90 

492 

44/5 

26 

4 

1 

y(2) 

MW-27  (4”) 

9/2/09 

1200 

100 

120/6 

32 

4 

3 

y(i) 

MW- 15  (2”) 

3/16/10 

100 

91 

9/4.5 

30 

2 

N/A 

y(i) 

MW-27  (4”) 

3/16/10 

340 

98 

70/5 

30 

4 

0 

n 

MW- 16  (2”) 

3/17/10 

105 

122 

13/5 

30 

2 

N/A 

y(i) 

MW-28  (4”) 

3/17/10 

165 

226 

37/5 

33 

4 

0 

n 

MW-24  (4”) 

8/16/10 

370 

302 

112/4.7 

33 

4 

N/A 

n 

4.4.5  Pour  Testing 

A  single  outdoor  field  test  of  methods  that  are  commonly  used  to  fill  VO  A  vials  was  conducted 
to  assess  the  degree  to  which  field  sample  handling  affects  analytical  results.  Three  different  fill 
rates  and  three  different  fill  methods  were  used.  The  fill  rates  for  the  tests  were  50  ml  per 
minute,  250  ml  per  minute,  and  1  L  per  minute.  The  fill  methods  included  a  vertical  “top  pour” 
approach  where  water  sample  was  allowed  to  free-fall  from  2  cm  above  the  top  of  the  vial  to  the 
bottom  of  the  vial;  a  “side  pour”  method  where  the  vial  was  held  at  a  30  to  45  degree  angle  to 
allow  the  discharge  to  flow  smoothly  down  the  side  of  the  vial;  and  a  “bottom  fill”  method  where 
the  discharge  tube  was  placed  at  the  bottom  of  the  VGA  vial  and  then  slowly  retracted  while 
keeping  the  end  of  the  tubing  slightly  submerged  during  the  filling  process  (Figure  18).  There 
were  six  replicates  for  each  combination  of  fill  rate  and  fill  method.  Sampling  was  first 
conducted  at  50  mL/min;  then  at  250  mL/min;  and  last  at  1  L/min.  It  took  approximately  50 
seconds  to  fill  a  vial  at  the  lowest  flow  rate,  approximately  15  seconds  at  250-mL/min  rate,  and 
approximately  3  seconds  at  1  L/min.  For  each  flow  rate,  the  sampling  order  was  reversed  with 
each  set  of  samples  (i.e.  a  set  of  samples  consisted  of  a  top-pour  sample,  a  side-pour  sample,  and 
a  bottom-fill  sample). 

Analytes  in  the  test  solution  included  benzene;  chloroform;  1,1-dichloroethane  (IIDCA);  1,2- 
dichloroethane  (12DCA);  1,1-dichloroethene  (IIDCE);  cis-l,2-dichloroethene  (cDCE); 
tetrachloroethene  (PCE);  toluene;  1,1,1-trichloroethane  (lllTCA);  trichloroethene  (TCE); 
trichlorofluoromethane  (TCEMA);  vinyl  chloride  (VC);  and  1,4-dioxane.  Sample  handling. 
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processing,  and  analyses  were  condueted  as  deseribed  previously.  Additional  information  on  the 
methods  used  can  be  found  in  Parker  and  Britt  (2012). 


Figure  18.  Photos  of  top-pour,  side-pour,  and  bottom- fill  methods,  respectively. 

For  each  of  the  twelve  analytes  found  in  the  test  solution,  the  concentration  for  a  given  fill 
method  and  fill  rate  was  divided  by  the  mean  eoncentration  of  that  analyte.  The  normalized 
values  were  then  pooled  into  one  data  set.  This  provided  us  with  a  data  set  with  233  values  for 
each  pour  method  or  fill  method  (or  699  total  values),  SigmaStat  (3.1)  software  (by  Systat 
Software  Inc.)  was  used  for  the  following  statistical  analyses.  A  two-way  Analysis  of  Variance 
(ANOVA)  was  conducted  to  determine  if  either  the  pour  method  or  the  flow  rate  had  an  effeet  on 
analyte  concentrations  for  the  entire  suite  of  analytes.  The  Holm-Sidak  method  for  a  pair-wise 
multiple  comparison  test  (Systat  Software,  Inc.  2004)  was  used  to  determine  which  flow  rates 
and  fill  methods  differed  from  each  other. 
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5.0  Results  and  Discussion 

5,1  Field  Investigation 

Several  phases  of  the  field  investigation  yielded  results  that  satisfied  the  goals  of  the  physieal 
field  testing  program.  Initial  multilevel  sampling  provided  information  that  allowed  seleetion  of 
appropriate  loeations  for  test  well  installation  at  the  field  site.  New  wells  largely  funetioned  as 
antieipated  and  allowed  eolleetion  of  aetive  and  passive  sampling  data  to  inform  our  researeh 
questions. 

Results  from  three  primary  foeus  areas  for  the  field  investigation  are  presented  here;  1) 
multilevel  passive  sampling  and  sensing  for  stratifieation  analysis,  2)  purge  sampling  with 
multiple  sample  eolleetion  for  “purge  eurve”  analysis,  and  3)  aneillary  testing  for  greater 
understanding  sample  eolleetion  variability. 

5.1.1  Stratification  Testing  Results 

Results  from  stratification  testing  of  the  existing  wells  indicated  a  stronger  degree  of 
stratification  during  the  phase  of  sampling  where  baffie  mixing  inhibitors  were  deployed 
between  passive  samplers.  Table  5  illustrates  results  for  several  constituents  in  the  isolated  vs 
un-isolated  sampling.  Consistent  with  previous  work  by  the  investigators,  there  were  differences 
in  degrees  of  stratification  with  the  alternate  deployment  configurations.  Somewhat  surprisingly, 
in  some  cases,  the  stratification  effect  inverted  using  the  zone  separation  method  and  was 
chemical-specific.  However,  not  all  examples  showed  substantial  differentials. 

TABLE  5,  Comparison  of  Open  Well  vs.  Baffle  Separator  Samples 


Well/Chemical 

BAFFLE 

maximum 

difference 

ratio* 

NO  BAFFLE 
maximum 
difference 
ratio 

Stratification 

Measurement 

improvement** 

MW9  (2  zones) 

Benzene 

-1.6 

-1.1 

45% 

11DCA 

1.7 

-1.2 

reverse 

11DCE 

1.1 

-1.1 

Reverse 

cis-1,2DCE 

1.1 

1 

10% 

PCE 

-3.3 

1.7 

Reverse 

TCE 

-1.4 

1.1 

Reverse 

VC 

3.7 

-3.6 

reverse 

MW14  (3  zones) 

Benzene 

21 

14 

50% 

11DCA 

144 

36 

300% 

11DCE 

28 

7.6 

268% 

cis-1,2DCE 

21 

9.2 

128% 

PCE 

-6.0 

-4.8 

25% 
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Well/Chemical 

BAFFLE 

maximum 

difference 

ratio* 

NO  BAFFLE 
maximum 
difference 
ratio 

Stratification 

Measurement 

improvement** 

TCE 

-1.6 

-2.5 

-36% 

VC 

129 

10 

1190% 

MW  17  (4  zones) 

Benzene 

~ 

~ 

11DCA 

2.0 

1.2 

67% 

11DCE 

-1.3 

1.5 

reverse 

cis-1,2DCE 

1.2 

-1.1 

reverse 

PCE 

~ 

~ 

TCE 

-1.4 

-1.2 

17% 

VC 

~ 

— 

MW21  (3  zones) 

Benzene 

- 

~ 

11DCA 

3.8 

7.1 

-46% 

11DCE 

2.5 

1.7 

47% 

cis-1,2DCE 

3.3 

2.7 

22% 

PCE 

1.9 

-1.3 

reverse 

TCE 

1.4 

-1.3 

reverse 

VC 

6.5 

5.1 

27% 

MW24  (4  zones) 

Benzene 

- 

~ 

11DCA 

1.5 

1.3 

15% 

11DCE 

-2.0 

-1.1 

82% 

cis-1,2DCE 

1.2 

1.1 

9% 

PCE 

-1.4 

1.3 

reverse 

TCE 

-1.9 

-1.3 

46% 

VC 

~ 

~ 

*  difference  between  highest  and  iowest  concentration,  negative  denotes 
concentrations  are  higher  at  the  bottom  of  the  weii 


**  positive  percentage  is  greater  stratification  present  with  zone  isoiation; 
negative  percentage  is  iesser  stratification,  reverse  denotes  stratification 
direction  changed  with  and  without  zone  isoiation 


These  preliminary  stratification  data  were  used  to  select  locations  where  stratification  was 
relatively  strong  for  installation  of  new  wells.  Three  locations  were  selected.  Figure  19  shows 
the  locations  relative  to  the  existing  stratification  test  wells.  As  described  in  the  Materials  and 
Methods  section,  these  wells  were  constructed  to  allow  multiple  sampling  approaches  with  a 
variety  of  downhole  sensors  and  sampling  equipment. 
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5.1.2  Purge  Curve  Testing  Results 

Using  these  new  wells  for  flow-through  dilution  tests,  purge  sampling,  and  multilevel  passive 
sampler/sensor  deployment,  multiple  test  seenarios  were  completed.  These  scenarios  included 
chemical  sensor  data  collection  during  a  pre-sampling  deployment  period;  pumping  with 
multiple  periodic  chemical  sample  collection;  and  pre-  and  post-purge  downhole  sample  capture. 

Purge  testing  data  show  that  pump  location  relative  to  contaminant  inflow  position  had  a 
sometimes  dramatic  effect  on  the  shape  of  the  purge/contaminant  curve.  That  is,  the  closer  or 
farther  from  the  contaminant  inflow  position  the  pump  was  located,  the  more  dramatic  the 
divergence  from  the  starting  and  ending  inflow  concentrations.  Figure  20(A),  shows  TCE 
concentration  drop  by  about  a  factor  of  3  in  early  pumping,  with  original  concentrations 
returning  after  about  three  well  volumes  pumped.  In  contrast.  Figure  20(B)  shows  a  smaller  50% 
boost  in  concentrations  early  in  the  purge  cycle.  The  implication  of  these  early  purge  changes  is 
that  chemical  concentration  changes  may  be  a  source  of  sampling  variability  that  is  not  widely 
recognized.  Further,  only  when  field  parameter  measurements  also  show  divergences  would  the 
field  sampling  practitioner  know  that  such  divergences  are  occurring.  Additional  results  below 
illustrate  that  such  data  is  not  necessarily  indicative  of  such  chemical  concentration  divergences. 
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Figure  20  (A/B).  Examples  of  purge  curves  for  TCE  in  MW-29  and  1,1 -DCE  in  MW-27.  Early  time  divergence 
from  flow  weighted  average  concentrations  during  purging.  Light  blue  and  pink  highlighted  samples  are  Snap 
Samples  collected  in  situ  above  (pink)  and  below  (blue)  the  pump  intake. 

Figure  21(A)  illustrates  a  ease  where  very  little  ehange  oeeurred  in  TCE  eoneentration 
throughout  the  purge,  until  the  very  end  when  the  purge  flow  rate  was  ehanged  (by  testing 
design).  Figure  21(B)  also  shows  minor  slow  ehange  in  eis-l,2-DCE  eoneentration  over  the 
eourse  of  the  purge,  but  a  fairly  rapid  ehange  in  eoneentration  when  the  purge  rate  was  ehanged 
at  the  end  of  the  purge  eyele.  Additionally,  the  Snap  Samples  above  and  below  the  pump  intake 
colleeted  at  the  beginning  of  the  purge  differ  substantially  from  those  concentrations  observed 
after  purging.  Despite  minimal  changes  in  VOC  concentration  during  the  course  of  the  purge, 
there  are  clearly  different  concentrations  of  VOC  contaminants  entering  the  well  from  different 
positions.  After  purging,  the  water  collected  with  the  Snap  Samplers  show  undetectable  VOC 
concentrations  above  the  pump  and  approximately  25%  higher  than  the  pump  discharge 
concentrations  below  the  pump  intake.  This  differential  can  be  interpreted  to  mean  that  the 
pump  discharge  contains  approximately  80%  VOC-containing  water  entering  the  well  from 
below  the  pump  position  and  20%  uncontaminated  water  entering  from  above  the  pump. 
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Figure  21  (A/B),  Examples  of  purge  curves  for  TCE  cis-l,2-DCE  in  MW -24.  Relatively  stable  concentrations 
throughout  purging  disrupted  by  changing  the  purge  rate  at  the  very  end  of  purging.  Light  blue  and  pink  highlighted 
samples  are  Snap  Samples  collected  in  situ  above  (pink)  and  below  (blue)  the  pump  intake. 

Appendix  A  contains  many  additional  data  sheets  from  the  several  “purge  curve”  tests  conducted 
during  the  investigation.  Additional  data  plots  are  presented  in  the  appendix  for  multiple 
chemicals  provide  additional  illustration  of  the  differences  among  test  scenarios  and  different 
wells. 

Another  important  observation  from  the  late  purge  data  comes  from  the  VOC  data  collected  in 
close  sequence  immediately  after  lowering  the  purge  rate.  Reducing  purge  rate  is  commonly 
recommended  to  improve  VOC  recovery  during  sample  collection  (USEPA  1985,  USEPA  2002, 
ASTM  2002).  When  contaminants  entering  the  well  bore  are  stratified,  changing  the  purge  rate 
changes  the  weighting  of  water  entering  the  pump  intake  because  water  level  changes.  When 
reducing  flow  rate,  water  level  rises  in  the  well.  This  water  level  change  causes  overweighting 
of  water  entering  below  the  pump.  In  the  case  shown  in  Figure  21,  this  results  in  a  temporary 
increase  in  contaminant  concentration. 

Lab  experiments  illustrate  a  similar  effect  with  the  opposite  condition — ^with  “contaminated” 
water  above  the  pump  intake.  Figures  22  and  23  show  the  photo  series  of  samples  collected 
before  and  after  changing  purge  rate  in  a  pipe  test.  Importantly,  the  change  in  purge  rate  can 
either  increase  or  decrease  pump  discharge  concentrations,  sometimes  dramatically.  In  the 
laboratory  case,  the  “contaminant”  entering  above  the  discharge  point  nearly  disappeared  for 
several  minutes  until  water  level  stabilized  and  the  flow-weighted  condition  returned.  For  field 
sampling,  this  effect  will  likely  go  unnoticed  and  could  be  a  substantial  source  of  sampling 
variability  if  sampling  is  conducted  immediately  after  lowering  the  flow  rate  (the  normal 
condition). 
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Drawdown  Stable  water  level,  stable  concentration 
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stable  pump  rate  down 

Figure  22.  Simulated  well  with  dyed  water  entering  above  pump  intake,  clear  water  entering  below  (out  of  view). 
Constant  head  reservoir  tank  with  clear  and  dyed  water  controls  equal  flow  into  upper  and  lower  parts  of  pipe. 
Discharge  in  controlled  by  gravity  with  a  valve  and  “pump”  position  located  equidistant  from  the  dyed  and  clear 
water  entry  points.  Flow  rate  is  1000ml  per  minute. 
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Figure  23.  Flow  rate  reduced  to  250ml  per  minute,  causing  rising  water  in  simulated  well.  Rising  water  causes 
discharge  to  temporarily  overweight  water  entering  from  the  bottom  half  of  the  well. 

Overall,  these  data  show  that  pump  position  can  be  an  important  factor  when  it  comes  to  rapidity 
of  attainment  of  flow-weighted  average  concentrations  during  purging,  and  the  degree  to  which 
the  discharge  changes  concentrations  throughout  purging.  Additionally,  the  data  suggest  that 
there  may  be  times  when  “purge  parameters”  such  as  temperature,  conductivity  or  oxidation 
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reduction  potential  may  mislead  a  sampling  team.  As  shown  in  Figure  24,  purge  parameter 
stabilization  may  appear  to  occur  at  a  time  when  contaminant  concentrations  have  not  yet 
stabilized  (parameter  data  from  same  purge  event  for  MW-24  as  Figure  21A/B  and  for  MW29  as 
Figure  20A).  Indeed,  apparent  stability  when  observing  reading-to-reading  results  at  a  purge 
volume  of  approximately  one  well  volume  pumped  may  actually  lead  the  sampling  team  to  collect 
samples  at  a  time  when  contaminant  concentrations  at  the  furthest  point  from  the  flow-weighted 
average  concentration  during  the  course  of  the  purge. 
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Figure  24,  Purge  curve  data  illustrate  small  and  slow  changes  in  parameter  readings. 

Alternately,  for  some  analytes,  stability  of  contaminant  concentrations  was  never  achieved, 
despite  long  completed  parameter  stability.  The  constituents  include  those  that  readily  degrade 
in  the  presence  of  oxygen  (e.g.  benzene,  acetone).  Figure  25  illustrates  examples  of  these 
compounds  in  the  MW29  purge  curve.  These  findings  may  indicate  the  well  itself  is  a  locus  of 
biodegradation.  In  these  circumstances,  passive,  sensor-based,  and  even  low-fiow/low-volume 
purge  approaches  may  be  problematic  for  representing  formation  groundwater  without  a 
thorough  purge.  A  technical  paper  is  in  preparation  that  includes  more  detail  on  results  and 
conclusions  of  this  aspect  of  the  study  (Britt  et  al,  in  prep).  These  findings  are  not  repeated  in 
detail  here  to  avoid  rejection  of  the  manuscript  as  “already  published”  material.  Appendix  A 
contains  additional  data  plots  of  the  purge/contaminant  curves.  Raw  data  lab  reports  are 
included  in  Appendix  C. 
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Figure  25.  Continued  increasing  concentrations  of  some  compounds  through  5  well  volumes  purged. 


The  tube  wells  installed  in  MW-27,  MW-28,  and  MW-29  (illustrated  in  Figures  15  and  16)  were 
utilized  to  assess  contaminant  concentrations  entering  the  well  but  not  yet  impacted  by  mixing 
and  flow  processes  within  the  well.  This  approach  was  somewhat  imperfect  in  the  ambient  flow 
regime  because  installation  difficulties  and  possible  flow  gradient  changes  during  the  course  of 
the  study  made  it  impossible  to  know  for  sure  if  the  tube  wells  were  always  on  the  upgradient 
(unaffected)  side  of  the  well.  Samples  collected  during  purging  were  less  likely  to  be  affected  by 
such  orientation  mismatch.  Some  of  the  purge  curve  analyses  included  collection  of  tube  well 
samples  at  the  end  of  the  large  volume  removal  actions.  Similar  to  those  interpretations 
allowable  by  collecting  samples  from  above  and  below  the  pump  intake,  tube  well  samples 
provided  additional  information  about  what  contaminants  were  entering  the  well  at  what 
position.  Figure  26  illustrates  the  range  of  tube  well  samples  compared  to  the  in-well  Snap 
Sampler  results  collected  from  above  and  below  the  pump  intake,  and  pumped  sample  results. 
Appendix  A  contains  many  more  examples  of  different  wells  and  different  contaminants. 
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Figure  26.  Purge  sampling  results  compared  to  Snap  Samples  collected  above  and  below  pump  intake  (pink  and 
medium  blue)  in  contract  to  tube  well  samples  collected  just  outside  the  well  (light  blue,  orange  and  purple 
highlights) 
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Tube  wells  were  also  used  in  a  passive  configuration  with  multiple-interval  passive  samplers.  To 
limit  the  effects  of  in-well  mixing  on  both  the  in-well  passive  devices,  but  also  the  indirect  effect 
on  the  tube  well  samples,  plastic  baffle  mixing  inhibitor  devices  were  placed  between  the  passive 
samplers  (Figure  13).  Polyethylene  Diffusion  Bag  (PDB)  samplers  were  deployed  at  1  foot 
intervals  in  two  of  the  wells  with  tube  wells.  Tube  well  samples  were  collected  before 
deployment  of  the  PDBs  and  then  just  prior  to  retrieval  of  the  PDBs.  There  are  only  4  tube  well 
samples,  while  there  were  10  PDB  samples  in  each  well.  The  tube  well  samples  reasonably 
matched  the  PDB  samples,  but  there  were  some  discrepancies.  Figure  27  illustrates  TCE  data  for 
MW-27  during  this  test.  Many  other  examples  are  included  in  Appendix  A,  Tables  A- 14  and  A- 
15. 


Figure  27.  Tube  well  samples  (Event  1)  collected  prior  to  insertion  of  the  baffled  PDB  samples;  tube  well  samples 
collected  just  prior  to  PDB  removal;  and  PDB  samples  collected  at  1  foot  intervals  within  the  screen  zone.  Stronger 
stratification  shown  in  the  Event  2  tube  wells  with  zone  isolation  present  in  the  well. 


5.1.3  Sample  Pour  Testing  Results 

In  addition  to  the  passive  sampling  and  purge  curve  sampling  analysis,  test  of  sample  pour 
methods  was  conducted  to  assess  the  degree  to  which  field  sample  handling  affects  analytical 
results.  Pour  tests  included  three  different  bottle  fill  rates  and  three  fill  methods.  Fill  rates  for 
the  tests  were  50ml  per  minute,  250ml  per  minute  and  1000  ml  per  minute.  Fill  methods 
included  a  vertical  “top  pour”  approach  where  water  was  introduced  to  the  VOA  vial  directly  by 
free-fall  from  2cm  above  top  of  the  vial  to  the  bottom;  a  “side  pour”  method  where  water  flow 
was  directed  to  the  side  of  the  bottle;  and  a  “bottom  fill”  method  where  the  discharge  tube  was 
placed  to  the  bottom  of  the  VOA  vial  and  kept  slightly  submerged  during  bottle  filling.  Six 
replicates  of  each  fill-rate  and  fill-method  combination  were  collected  and  analyzed  for  VOCs. 

The  results  for  each  of  the  pour  test  methods  for  the  50-mL/min,  250-mL/min,  and  1-L/min  flow 
rates  are  summarized  in  Table  6.  The  raw  data  can  be  found  in  Appendix  A  Table  A-15. 
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Table  6.  Summary  of  the  pour-test  study  results. 


Analvte 

50BF 

50SP 

50TP 

250BF 

250SP 

250TP 

lOOOBF 

lOOOSP 

lOOOTP 

Benzene 

Mean 

Cone.* 

0.88 

0.87 

0.81 

0.95 

0.88 

0.87 

0.91 

0.89 

0.88 

Std.  dev. 

0.069 

0.029 

0.045 

0.074 

0.049 

0.057 

0.024 

0.19 

0.13 

%  RSD 

7.8 

3.4 

5.6 

7.8 

5.6 

6.6 

2.7 

21.7 

14.3 

Chloroform 

Mean 

Cone.* 

21.7 

20 

18.6 

22.7 

20.3 

20.4 

22.9 

21.5 

22.1 

Std.  dev. 

1.28 

0.98 

0.61 

0.82 

4.34 

0.71 

0.48 

3.6 

2.02 

%  RSD 

5.9 

4.9 

3.3 

3.6 

21.4 

3.5 

2.1 

16.7 

9.1 

IIDCA 

Mean 

Cone.* 

99.9 

92 

85.1 

105 

93 

91.6 

106 

97.6 

100 

Std.  dev. 

5.74 

5.36 

2.91 

5.91 

20.6 

5.7 

3.4 

19.1 

10.7 

%  RSD 

5.7 

5.8 

3.4 

5.6 

22.2 

6.2 

3.3 

19.6 

10.7 

12DCA 

Mean 

Cone.* 

2.35 

2.27 

2.18 

2.43 

2.34 

2.34 

2.58 

2.54 

2.58 

Std.  dev. 

0.151 

0.063 

0.117 

0.39 

0.253 

0.136 

0.101 

0.292 

0.121 

%  RSD 

6.4 

2.8 

5.3 

16.1 

10.8 

5.8 

3.9 

11.5 

4.7 

Toluene 

Mean 

Cone.* 

0.83 

0.81 

0.76 

0.88 

0.83 

0.84 

0.87 

0.82 

0.84 

Std.  dev. 

0.026 

0.045 

0.014 

0.036 

0.043 

0.022 

0.039 

0.123 

0.077 

%  RSD 

3.2 

5.6 

1.9 

4.1 

5.2 

2.6 

4.6 

15 

9.2 

lllTCA 

Mean 

Cone.* 

1.62 

1.36 

1.44 

1.54 

1.75 

1.36 

1.79 

1.61 

1.67 

Std.  dev. 

0.104 

0.38 

0.085 

0.442 

0.102 

0.462 

0.126 

0.481 

0.261 

%  RSD 

6.4 

27.9 

5.9 

28.6 

5.8 

33.9 

7 

29.9 

15.7 

TCE 

Mean 

Cone.* 

120 

112 

104 

123 

105 

110 

121 

113 

114 

Std.  dev. 

2.61 

2.71 

2.32 

2.8 

30.3 

1.79 

3 

26.2 

13.5 

%  RSD 

2.2 

2.4 

2.2 

2.3 

28.8 

1.6 

2.5 

23.2 

11.9 

TCFMA 

Mean 

Cone.* 

48.2 

44.1 

40.2 

49 

40.5 

40.7 

49.7 

42.8 

44.3 

Std.  dev. 

3.36 

2.8 

1.37 

5.25 

16.75 

5.82 

2.35 

15.15 

9.71 

%  RSD 

7 

6.4 

3.4 

10.7 

41.4 

14.3 

4.7 

35.4 

21.9 

IIDCE 

Mean 

Cone.* 

239 

218 

198 

245 

206 

207 

242 

214 

221 

Std.  dev. 

17.3 

16.9 

6.9 

22.1 

66.9 

10.1 

11.1 

63.4 

41.1 

%  RSD 

7.2 

7.8 

3.5 

9 

32.5 

4.9 

4.6 

29.7 

18.6 

cDCE 

Mean 

Cone.* 

49.7 

45.7 

42.4 

51.3 

45.8 

45.5 

51.2 

48.1 

49.2 

Std.  dev. 

3.22 

2.83 

1.54 

2.48 

8.7 

1.71 

1.42 

8.27 

4.34 

%  RSD 

6.5 

6.2 

3.6 

4.8 

19 

3.8 

2.8 

17.2 

8.8 

tDCE 

Mean 

Cone.* 

BDL** 

BDL** 

BDL** 

BDL** 

BDL** 

BDL** 

0.72 

0.74 

0.76 

Std.  dev. 

0.202 

0.206 

0.155 

%  RSD 

28.1 

27.9 

20.3 

PCE 

Mean 

Cone.* 

130 

123 

114 

132 

111 

119 

129 

117 

119 

Std.  dev. 

2.28 

3.35 

1.86 

2.1 

37.5 

3.1 

2.1 

29.8 

16.8 

%  RSD 

1.8 

2.7 

1.6 

1.6 

33.9 

2.6 

1.6 

25.4 

14.2 
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Analvte 

50BF 

50SP 

50TP 

250BF 

250SP 

250TP 

lOOOBF 

lOOOSP 

lOOOTP 

Mean 

Cone.* 

26.1 

25 

22.4 

30 

23.8 

25.4 

29.4 

26.7 

26.4 

Std.  dev. 

3.94 

2.07 

0.84 

2.65 

6.64 

1.58 

0.88 

8.53 

6.28 

%  RSD 

15.1 

8.3 

3.8 

00 

bo 

28 

6.2 

3 

31.9 

23.8 

1  An-  1  Mean 

14Dioxane  „ 

Cone.* 

1122 

1186 

1138 

1172 

1174 

till 

1317 

1310 

1367 

Std.  dev. 

89.1 

154.7 

100.1 

106 

124 

96 

145 

43.6 

45.1 

%  RSD 

7.9 

13 

8.8 

9.1 

10.6 

8.7 

11 

3.3 

3.3 

*  Concentration  in  |tg/L 

**  BDL,  most  below  detection  level 

’  Minus  one  outlier 

The  statistical  analyses  showed  that  the  differences  among  the  different  fill  rates  and  bottle-fill 
methods  were  highly  significant  and  that  the  interaction  of  the  fill  rate  and  fill  method  was  also 
significant.  This  significant  interaction  means  that  the  effect  of  fill  rate  differed  among  the  fill 
methods,  and  vice  versa  (i.e.,  the  effect  of  fill  method  differed  at  different  fill  rates).  Therefore, 
the  impact  of  fill  rate  and  fill  technique  were  examined  separately. 

5. 1.3.1  Fill  Method 

There  was  a  significant  difference  between  the  bottom-fill  and  side-  pour  methods  and  between 
the  bottom-fill  method  and  the  top-pour  method  at  all  three  fill  rates  (Table  7).  For  these 
comparisons,  a  number  followed  by  two  letters  were  used  to  differentiate  among  the  methods. 
The  number  gives  the  vial  fill  rate  in  milliliters  per  minute  and  the  letters  describe  the  fill 
method.  For  example,  50BF  means  that  the  vial  was  filled  at  50  mL/min  using  the  bottom-fill 
method.  The  number  1000  is  used  to  denote  the  fastest  flow  rate,  -1000  mL/min  (-1  L/min.). 
The  “>”  symbol  means  that  the  difference  in  VOC  recovery  was  greater  at  statistically  significant 
level,  and  the  symbol  representing  no  significant  difference  in  recovery. 
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Table  7,  Summary  of  the  statistical  analyses  comparing  the  different 
fill  methods 


%  Difference 

Significant? 

P= 

Fill  Method  (BF  vs.  SP) 

50BF  >  50SP 

7.5 

Yes 

0.005 

250BF  >  250SP 

14.3 

Yes 

<0.001 

1000BF  >  1000SP 

8.1 

Yes 

0.021 

Fill  Method  (SP  vs.  TP) 

50SP  >  50TP 

6.1 

Yes 

0.021 

250SP  ~  250TP 

3.3 

No 

0.199 

1000SP~  1000TP 

2.3 

No 

0.358 

Fill  Method  (BF  vs.  TP) 

50BF  >  50TP 

13.6 

Yes 

<0.001 

250BF  >  250TP 

11.4 

Yes 

<0.001 

1000BF  >  1000TP 

5.7 

Yes 

0.001 

In  both  cases,  recovery  was  best  (i.e.,  concentrations  were  highest)  with  the  bottom-fill  method, 
and  recovery  was  approximately  6  to  14%  greater  than  with  either  the  side -pour  method  or  the 
top-pour  method  at  all  three  fdl  rates.  When  the  side-pour  and  top-pour  methods  were 
compared,  there  was  no  signifieant  differenee  between  the  two  methods  at  the  two  faster  flow 
rates  but  there  was  at  the  slowest  fill  rate. 

The  rankings  for  the  fill  methods  for  each  of  the  flow  rates  are  summarized  here; 

50BF  >  50SP  >  50TP 
250BF>250SP~250TP 
lOOOBF  >  1000  SP  ~  lOOOTP 

These  results  are  contrary  to  guidance  documents  that  recommend  a  side-pour  method.  Also, 
contrary  to  conventional  practice,  this  study  shows  that  the  side-pour  method  is  no  better  than 
simply  pouring  straight  down  into  the  vial  from  the  top,  except  when  a  very  low  fill  rate  is  used. 

5. 1.3. 2  Fill  Rate 

There  were  significant  differences  in  the  recoveries  of  the  VOCs  between  the  fastest  and  slowest 
fill  rates  for  all  three  fill  methods  (Table  8).  In  all  cases,  recovery  was  greater  (i.e., 
concentrations  were  higher)  with  the  faster  flow  rate,  and  the  differences  in  recovery  ranged 
from  ~6  to  14%  for  the  three  fill  methods.  When  the  250-mL/min  and  1-L/min  flow  rates  were 
compared,  there  was  a  significant  difference  in  the  reeoveries  for  the  side-pour  method  and  for 
the  top-pour  method,  with  recovery  7  to  8%  greater  for  the  faster  flow  rate.  In  contrast,  when  the 
recovery  of  the  higher  two  flow  rates  was  compared  for  the  bottom-fill  method,  the  pereent 
difference  in  recovery  was  only  1 .6%  and  this  difference  was  not  statistieally  significant.  When 
the  two  lower  flow  rates  were  eompared,  there  was  no  significant  difference  in  the  reeovery  for 
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the  bottom-fill  method  or  for  the  side-pour  method  but  there  was  for  the  top-pour  method.  For 
the  top-pour  method,  recovery  was  7%  greater  at  250  mL/min  than  at  50  mL/min. 

Table  8,  Summary  of  the  statistical  analyses  comparing  the  different  fill  rates. 


Fill  Rate 

50  VS.250  mL/min 

%  Difference 

Significant? 

P= 

50BF  ~  250BF 

4.9 

No 

0.04 

50SP  ~  250SP 

2.2 

No 

0.404 

50TP  <  250TP 

7.1 

Yes 

0.006 

250  vs. 1000  mL/min 

250BF-1000BF 

1.6 

No 

0.499 

250SP  < 1000SP 

8.2 

Yes 

0.002 

250TP  < 1000TP 

7.3 

Yes 

0.006 

50  vs.  1000  mL/min 

50BF  <  1000BF 

6.5 

Yes 

0.007 

50SP  <  1000SP 

6 

Yes 

0.023 

50TP  <  1000TP 

14.6 

Yes 

<0.001 

The  rankings  for  the  recoveries  for  each  of  the  pour  methods  at  each  of  the  flow  rates  are 
summarized  below. 

1000  BF~  250  BF;  250BF  ~  50BF;  lOOOBF  >  50BF 

1000SP>250SP~50SP 

1000TP>250TP>50TP 

These  results  are  contrary  to  conventional  wisdom  and  guidance  that  recommends  using  slower 
fill  rates  (e.g.,  100  mL/min  or  less)  to  reduce  turbulence  (USEPA  2002,  ASTM  2002).  Also, 
contrary  to  conventional  wisdom,  turbulence  due  to  fast  vial  filling  did  not  appear  to  be  a  factor 
in  VOC  losses  as  vigorous  mixing  was  evident. 

Sampling  order  with  respect  to  the  fill  rate  used  was  not  randomized  in  the  test  (i.e.,  the  bottles 
were  filled  first  at  50-mL/min,  then  at  ~  250-mL/min,  and  finally  at  --l-L/min),  therefore  it  is 
possible  that  this  affected  the  recovery  of  the  analytes  in  some  fashion.  Flowever,  sampling  just 
a  few  liters  from  a  fixed  position  in  a  55  gallon  drum,  a  mechanism  for  time  variation  of 
concentration  during  sampling  is  not  evident,  but  can’t  be  completely  ruled  out.  Parker  and  Britt 
(2012)  contains  additional  details  and  explanation  of  these  results  and  concerns. 

5. 1.3. 3  Variability 

There  was  also  variability  in  the  data  for  the  different  fill  methods  and  flow  rates.  This 
information  is  summarized  in  Table  9,  which  gives  the  pooled  Percent  Relative  Standard 
Deviation  for  each  of  the  sampling  methods  and  flow  rates.  The  variability  was  generally  less  at 
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the  lower  fill  rate.  This  is  especially  true  for  the  side -pour  and  top-pour  methods.  Overall,  the 
bottom-fill  method  had  the  least  variability.  In  contrast,  variability  was  greatest  when  the  side- 
pour  method  was  used,  especially  at  the  two  higher  flow  rates.  This  finding  also  suggests  the 
Bottom  Fill  method  is  an  improvement  over  standard  side  pour  approach. 


Table  9,  Pooled  %RSD  values  for  each  fill  method  and  flow  rate. 


Pooled  Percent  Relative  Standard  Deviation 

Flow  rate 

BF 

SP 

TP 

1  L/min 

2.14 

10.4 

6.3 

250  mUmin 

4.94 

10.5 

5.04 

50  mUmin 

3.23 

4.43 

2.05 

5. 1.3. 4  Discussion 

Contrary  to  conventional  thought,  the  results  of  this  analysis  suggest  that  faster  fill  rates 
generally  yield  higher  VOC  recoveries.  Also  contrary  to  conventional  thought  on  the  subject, 
pouring  down  the  side  of  the  sample  container  is  not  substantially  better  than  simply  pouring 
down  from  the  top  of  the  vial  with  sample  in  free  fall.  Clearly  the  best  pour  method  at  the  lower 
fill  rates  was  the  bottom  fill  approach.  At  high  fill  rates  the  fill  method  was  less  important, 
apparently  because  the  sample  vial  was  filled  in  just  a  few  seconds  and  exposure  was  minimized. 
The  results  of  this  test  strongly  suggest  that  “turbulence”  of  the  sample  water  was  not  a  driver  for 
VOC  loss,  but  rather  time  of  exposure  and  surface  area  of  exposure  drive  VOC  loss.  Slow  bottle 
fill  rates  such  as  50ml  per  minute  and  pouring  down  the  side  of  the  vial  or  from  the  top  yielded 
the  lowest  VOC  recoveries.  Bottom  filling  at  low  flow  rate  significant  improved  recovery. 
Bottom  filling  at  the  250ml  flow  rate  also  significantly  improved  recovery  over  top  pour  and  side 
pour  methods.  Based  on  these  results,  the  bottom  fill  approach  should  normally  be  advised  at 
low  to  moderate  fill  rates,  and  higher  fill  rates  should  be  employed  when  possible  and  within 
reason.  Physical  constraints  exist  for  VGA  vial  filling  that  may  make  rates  approaching  and 
above  lOOOml/min  more  difficult.  However,  our  testing  showed  that  fill  rates  up  to  lOOOml/min 
using  the  bottom  fill  approach  were  achievable.  The  key  point  being  that  reduction  of  flow  rate 
to  accommodate  vial  filling  is  not  needed  at  rates  where  the  practicality  of  vial  filling  is  not 
problematic  (i.e.  there  is  no  reason  to  reduce  flow  from  250ml/minute  to  50  or  lOOml/min  to 
collect  samples). 


2014 
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5,2  Analytical  and  Numerical  Models  of  Monitoring  Well  Pumping 

Many  of  the  purging  curves  in  the  field  tests  showed  deviations  similar  to  those  predicted  by  the 
analytical  models  of  Martin-Hayden  (2000a).  Those  analytical  models  were  based  on  several 
simplifying  assumptions  including,  inviscid  flow,  predominantly  vertical  flow  within  the  well 
bore,  and  pumps  positioned  at  either  end  of  the  screen.  In  order  to  account  for  the  effects  of 
viscosity,  radial  flow  and  variable  pump  positions,  more  complex  analytical  and  numerical 
models  were  derived  for  this  project.  Analytical  and  numerical  models  were  developed  to 
account  for  Poisueille  (viscous  pipe)  flow  and  investigate  the  influences  of  viscosity  on  flow 
paths,  travel  time  and  concentration  distributions. 

The  models  indicate  that  groundwater  entering  along  the  well  screen  during  pumping  (Figure 
28a)  will  accumulate  and  velocities  will  accelerate  as  the  well  water  approaches  the  pump  intake 
(Figure  28c).  The  flow  paths  and  travel  times  that  result  from  this  velocity  distribution  are  used 
to  calculate  (1)  the  advective  front  of  groundwater  entering  the  well  screen  as  a  function  of  time, 
(2)  the  time  it  takes  groundwater  to  travel  from  the  screen  to  the  pump  intake,  (3)  the  zone 
captured  by  the  pump  intake  giving  the  position  of  origin  of  the  water  entering  the  pump  intake, 
and  (4)  pumping  curves  giving  the  fraction  of  pumped  water  that  is  initial  well  water  as  pumping 


(a)  (b)  (c)  (d) 

Figure  28.  (a)  Groundwater  entering  along  screen,  (b)  Parabolic  “Poiseuille”  velocity  distribution  and  stream 
lines,  (c)  Accelerating  velocities  approaching  pump  intake,  (d)  advancing  advective  front  or  annulus  of 
groundwater  (in  violet) 
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proceeds.  Here  the  volume  of  well  water  is  eonsidered  the  volume  of  water  within  the  sereen 
seetion  of  the  well  (Vw)  that  is  being  displaced  by  inflowing  groundwater  and  volume  pumped 
(Vp)  is  calculated  relative  to  that  volume,  i.e.,  Nwv=VpA^w 

(1)  The  advective  front  is  essentially  an  annulus 

of  in-flowing  groundwater  that  is  displacing 

pre-pumping  well  water  of  radius  r  (shown 

on  Fig.  17d  by  the  fine  white  dashed  lines 

and  violet  zone)  that  advances  at  a  rate  of 
1/2 

r=[(l-tanh(Nwv)]’  •  Because  the  adveetive 
front  is  independent  of  the  volume  pumped, 
the  remaining  well  water  is  a  eylinder  of 
uniform  radius  r.  This  eharaeteristic  of  well 
water  displacement  during  pumping  is  also 
demonstrated  by  physical  and  numerical 
models  of  monitoring  well  pumping  (shown 
as  the  warm  colors  on  Figure  29a  and  29b). 

Here  the  well  is  pumped  at  the  bottom  and 
the  annulus  of  groundwater  shown  as  blue  is 
replaeing  the  cylinder  of  well  water  shown 
as  red.  Mixing  is  indicated  by  the  diffuse 
zone  between  the  well  water  and 
groundwater. 

(2)  The  travel-time  between  a  point  along  the 
screen  and  the  pump  intake  is  given  by 
Nwv=cosh'  (zo'  )  where  zi  is  the  entry 
position  (e.g.,  labels  a  through  f  on  Figure 
28).  Groundwater  entering  at 
position  b  along  the  screen  will 
reach  the  pump  intake  after  one 
quarter  of  a  well  volume  is 
pumped  (Nwv=0.25)  whereas 
groundwater  entering  at  point  d 
will  arrive  at  Nwv=l  • 

(3)  The  travel  time  function  is  also 
used  to  calculate  how  the  zone 
of  groundwater  and  well  water 
captured  by  the  pump  intake 
propagates  down  the  well,  i.e., 
z=sech  (Nwv).  For  example,  if 
half  of  a  well  volume  is  pumped 


Physical  Model  Numerical  Model 


(a)  (b) 


Figure  29.  (a)  Physical  and  (b)  numerical  models  of 
the  diffuse  annulus  of  groundwater  resulting  from 
uniform  inflow  and  mixing  at  the  screen. 
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the  groundwater  and  well  water  are  originated  from  the  level  that  eorresponds  to  (b).  If  the 
well  water  and/or  the  groundwater  are  heterogeneous,  eoneentrations  arriving  at  the  pump 
intake  will  fluetuate  as  a  funetion  of  this  travel  time  (e.g.,  see  the  diseussion  related  to  Figure 
30). 

(4)  The  volume  of  pre-pumping  well  water  remaining,  Vr  (relative  to  the  initial  volume  of  well 
water  Vw)  after  a  certain  volume  of  water,  N^v,  is  pumped  equals  the  relative  contribution  of 
well  water  at  the  pump  intake.  Fig.  M3  shows  various  theoretical  curves  of  pre-pumping 
well  water  remaining  as  a  function  of  volume  pumped  Nwv.  The  curve  labeled  “(F)  Exp(- 
Nwv)”  is  calculated  under  the  unrealistic  assumption  of  inviscid  (a.k.a.,  Darcian)  flow  within 
the  well.  This  corresponds  to  the  numerical  model  of  inviscid  flow  shown  as  the  curve 
labeled  “(A)  Inviscid”  and,  interestingly,  purging  a  well  with  thorough  mixing.  Because  a 
well  is  an  open  column  of  water  rather  than  a  porous  (Darcian)  medium,  the  velocity 
distribution  has  been  demonstrated  to  be  “Poiseuillian”  (parabolic  in  lateral  cross  section,  see 
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Figure  31.  (a)  concentration  distribution  within  the  monitoring  well  before  pumping  given  various  degrees  of 
mixing  as  groundwater  flows  through  the  well;  Blue:  no  mixing  of  through-flowing  groundwater,  i.e.,  well  water 
concentration  is  the  same  as  the  groundwater  concentration  distribution  adjacent  to  the  well,  Green:  partial  mixing. 
Orange:  thorough  mixing,  (b)  flow  regime  and  travel-time  model,  (c)  pumping  curves  for  various  pump  intake 
positions  (Pz)  relative  to  the  concentration  distribution,  i.e.  for  Pz=0  the  pump  intake  is  at  the  high  concentration 
end  of  the  screen,  Pz=50  at  the  midpoint  and  Pz=100  at  the  low  concentration  end.  Each  of  three  curves  in  a  set  of 
the  same  color  represents  one  of  the  pre-pumping  mixing  conditions. 
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Figure  28b).  Including  Poiseuille  flow  in  the  analytical  model  gives  C=sech  (N^v)  and  the 

curve  labeled  “(C)  Poiseuille.” 

In  the  presence  of  a  concentration  gradient  adjacent  to  the  well  screen  (Figure  31a)  the 
concentrations  arriving  at  the  pump  intake  will  vary  as  a  function  of  volume  pumped  (N^v)-  hi 
Figure  31a  the  groundwater  concentration  is  high  at  the  top  of  the  well  (blue  curve  on  Figure 
31a).  If  the  well  water  is  thoroughly  mixed  before  pumping  begins,  and  the  well  is  then  pumped 
from  the  top  (dark  blue  solid  line,  Pz=0,  Co=mixed)  then  the  concentrations  will  rise  rapidly 
early  (Nwv<0.25),  appear  to  stabilize  before  Nwv=l  and  gradually  approach  the  average 
groundwater  concentration.  If  groundwater  is  flowing  horizontally  through  the  well  and  remains 
unmixed  before  pumping  begins,  the  concentration  distribution  will  mimic  the  groundwater 
concentration  distribution  (i.e.,  stratification).  In  this  case,  during  early  pumping  the 
concentration  at  the  pump  intake  will  be  high  because  both  the  well  water  and  the  approaching 
groundwater  are  high  (dot-dashed  dark  blue  line  at  the  top  of  Figure  31c).  The  concentration 
will  then  decrease  rapidly  as  groundwater  and  well  water  of  lower  concentrations  arrive  at  the 
pump  intake.  If  the  well  water  is  partially  mixed,  the  well  water  concentration  will  be  irregular 
and  the  concentrations  arriving  at  the  pump  intake  will  fluctuate  between  the  unmixed  and 
thoroughly  mixed  condition  (the  irregular  dark  blue  line  on  Figure  31c).  These  three  degrees  of 
mixing  are  simulated  for  each  level  of  pump  intake  down  to  Pz=100  (the  dark  green  curves  on 
Figure  31  representing  pumping  from  the  bottom  of  the  well. 

Effluent  concentration  during  pumping  is  also  affected  by  the  position  of  the  pump  intake 
relative  to  contaminant  sources  in  the  surrounding  formation,  and  by  conductivity  variations  in 
the  formation.  Numerical  modeling  provides  a  means  of  examining  how  such  variations  affect 
effluent  concentration  and  the  time  to  reach  a  flow-weighted  average  concentration.  Figure  32a 
illustrates  a  system  with  a  strong  variation  in  formation  concentrations,  with  highest 
concentrations  located  in  a  zone  of  lOx  higher  hydraulic  conductivity  located  in  the  upper  part  of 
the  column  (between  dark  horizontal  lines  at  z=85  to  95).  The  pre-pumping  well  water 
concentration  is  a  thoroughly  mixed,  flow-weighted  average  of  the  concentrations  adjacent  to  the 
well.  During  pumping  from  near  the  bottom  of  the  screen  (Figure  32,  Pz=8  where  the  stream 
lines  converge)  the  concentration  starts  at  the  flow-weighted  average  since  the  first  water  to 
arrive  is  the  thoroughly  mixed  pre-pumping  well  water.  As  with  the  analytical  models  (Figure 
31c),  the  concentrations  rapidly  diverge  from  the  initial  concentration  during  early  pumping 
(Figure  32b),  with  the  direction  of  change  dependent  on  the  intake  position  relative  to  the 
contamination  zone.  With  the  pump  intake  located  far  from  the  solute  source,  effluent 
concentrations  first  decrease  during  pumping,  as  clean  water  is  drawn  rapidly  to  the  intake  to 
replace  the  initially  well-mixed  water  assumed  to  exist  in  the  well.  Only  later,  when  solutes  have 
had  sufficient  time  to  advect  from  the  more  distant  solute  source,  does  the  concentration  revert  to 
the  flow- weighted  average. 
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This  effect  is  more  pronounced  in  the  presence  of  a  hydraulic  conductivity  heterogeneity.  Also 
similar  to  the  analytical  models,  the  concentration  gradually  approaches  the  flow  weighted 
average  concentration  at  late  times.  Major  benefits  of  numerical  models  include  the  ability  to 
include  more  complex  relationships  between  hydraulic  and  concentration  heterogeneities,  the 
influences  transient  effects  such  as  solutal  and  thermal  convection,  and  affects  of  changing 
pumping  conditions.  As  further  field  data  is  collected  on  these  relationships,  further 
complexities  will  be  included  in  the  numerical  modeling  to  quantify  the  effects  of  these 
processes. 


Figure  32.  Results  of  numerical  simulation  of  pumping  from  a  fully  screened  well  in  a  formation  of  heterogeneous 
concentration  and  permeability.  Panel  a  displays  streamlines  and  concentration  distribution  after  30  minutes  of 
pumping  at  250  ml/hour  with  the  pump  intake  at  the  bottom  of  the  well  and  solute  source  in  a  layer  near  the  top  of 
the  domain.  Panel  b  shows  effluent  concentration  vs.  time  resulting  from  pumping  at  the  bottom  of  the  well,  Pz=  8 
(blue  line),  in  the  middle  of  the  screen  section  (green  line)  and  near  the  top,  Pz=88  (red  line). 


5.2.1  Flow  distortions  during  ambient  conditions 

Groundwater  will  flow  horizontally  across  the  well  only  if  hydraulic  gradients  are  strictly 
horizontal,  and  no  density- driven  mixing  occurs.  Because  the  well  presents  a  highly  permeable 
conduit  relative  to  the  permeability  of  the  surrounding  formation,  slight  vertical  hydraulic 
gradients  will  induce  vertical  wellbore  flow.  In  the  example  shown  in  Figure  33  a  vertical 
gradient  that  is  10%  of  the  horizontal  gradient  drives  a  slight  upward  flow  in  the  formation. 
When  that  slight  upward  flow  encounters  the  highly  conductive  wellbore  the  flow  lines  will  be 
refracted  sharply  upward  (Figure  33  red,  yellow  and  black  arrows).  The  flow  lines  are  shown  by 
bands  of  different  colors  with  detail  at  the  top  and  bottom  of  the  screened  section  where  flow  is 
most  highly  distorted.  Inflow  is  focused  on  the  high-head  end  of  the  well,  the  bottom  of  the  well 
in  this  case,  and  attenuated  at  the  low-head  end. 
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The  focused  flow  at  the  high-head  end  of  the  screen 
will  be  spread  over  a  longer  segment  of  the  screen, 
whereas,  the  distributed  inflow  at  the  low  head  end 
will  be  focused  over  a  smaller  length  as  the  water 
leaves  the  well.  Because  groundwater  inflow  is 
focused  at  the  high  head  end  of  the  screen,  flow- 
weighted  average  concentrations  entering  the  well 
will  also  be  weighted  toward  that  end.  Within  the 
well  this  groundwater  will  be  transported  vertically 
and  most  likely  be  mixed  by  the  many  processes 
investigated  in  this  study.  Consequently  the 
concentrations  leaving  the  well  down  gradient  will 
altered  and  spread  over  a  larger  vertical  interval, 
i.e.,  the  shadow  zone. 

The  vertical  wellbore  flow  is  proportional  to  the 
vertical  gradient  in  the  surrounding  formation.  In 
the  limit,  as  the  vertical  component  of  the  gradient 
exceeds  the  horizontal  component,  vertical 
wellbore  flow  will  become  predominantly  vertical 
and  inflow  will  be  highly  concentrated  at  the  high- 
head  end  of  the  screen.  In  this  case  the 
groundwater  will  enter  the  high-head  end  of  the 
screen  and  exit  the  low-head  end.  The  well  water 
will  then  be  representative  of  the  zone  surrounding 
the  low-head  end  and  the  shadow  zone  will  be 
emanating  from  the  high-head  end. 

Because  vertical  gradients  tend  to  be  pronounced  in  recharge  and  discharge  areas  where  they 
tend  to  vary  in  time,  these  vertical  wellbore  flow  effects  are  expected  to  be  variable  in  intensity 
and  direction.  The  velocity  of  wellbore  flow  will  be  several  times  that  of  the  average  linear 
velocity  within  the  formation  due  to  the  focusing  of  flow.  For  example,  in  unconsolidated 
formations,  e.g.,  silt  and  sand,  where  typical  hydraulic  gradients  are  on  the  order  of  0.01  m/m, 
average  linear  velocities  would  be  between  millimeters  and  centimeters  per  day.  Thus,  wellbore 
flow  will  be  in  that  range  as  well,  significantly  slower  than  solutal  and  thermal  convection  and 
pumping  (even  for  pumping  during  low-flow  sampling,  e.g.,  1  L/min).  This  would  suggest  that 
groundwater  flowing  through  a  well,  regardless  of  the  vertical  or  horizontal  trajectory  within  the 
well,  would  likely  mix  in  most  cases. 
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5,3  Physical  Modeling 


5.3. 1  2D  and  3D  Physical  Models 


The  physical  model  of  monitoring  wells 
and  the  surrounding  environment  was 
constructed  to  be  flexible  in  terms  of 
horizontal  and  vertical  hydraulic  gradients, 
thermal  gradients  and  visibility  (Figure  34). 
The  model  is  essentially  a  rectangular  tank 
(length  x=90cm,  width  y=60cm,  height 
z=  120cm)  of  1  inch  thick,  transparent 
polycarbonate,  filled  with  fine  sand,  with 
specified  head  reservoirs  of  gravel  at  each 
end  (x=0  and  x=90cm).  The  hydraulic 
conductivity  of  the  gravel  (Kg=  1.2x10’' 
cm/s)  is  approximately  two  orders  of 
magnitude  greater  that  of  the  sand 
(Ks=3.8xl0’^  cm/s)  and  acts  to  vertically 
distribute  the  head  from  the  three  levels  of 


ports  on  each  reservoir.  By  adjusting  these 
heads,  a  vertical  gradient  is  imposed  across 
the  model.  The  slight  vertical  trajectory  of 
the  flow  lines  shown  on  Figure  34  would  Figure  34.  Physical  model  of  Half  Well  and  Full  3D  Well 
result  from  a  10%  vertical  gradient  with  controllable  horizontal  and  vertical  gradients. 

equivalent  to  those  used  to  produce  the  - 

flow  regime  in  Figure  33.  Adjusting  the  vertical  gradient  imposed  by  the  specified  head 
reservoirs  allows  the  investigation  of  vertical  flow  within  the  two  well  models,  i.e.,  the  partially 
penetrating  “Flalf  Well”  and  the  fully  penetrating  “Full  3D  Well”  shown  in  Figure  34. 


The  Half  Well  allows  the  flow  and  mixing  within  the  well  to  be  directly  observed  using  tracer 
dyes  and  recorded  with  time-lapse  photography,  e.g.,  the  method  used  to  produce  the  image  in 
Figure  29a  and  the  images  for  the  2-D  physical  models.  The  image  of  Figure  29a  was  produced 
by  equating  an  RGB  color  value  with  a  relative  concentration  value  stored  as  a  matrix.  This 
matrix  of  relative  concentrations  allows  the  mixing  to  be  characterized  by  comparing  the  results 
to  the  numerical  models,  e.g.,  the  comparison  of  Figure  29a  to  29b.  Figure  29a  was  produced  by 
introducing  a  uniform  distribution  of  red  dye  to  the  wellbore  before  pumping  and  taking  time- 
lapse  photos  as  pumping  proceeded.  As  the  un-dyed  water  entered  the  screen  it  produced  a 
diffuse  annulus  of  groundwater  as  predicted  by  the  analytical  models  and  corroborated  by  the 
numerical  models.  The  diffuse  nature  of  the  annulus  indicates  that,  as  shown  by  the  numerical 
models  of  dispersive  mixing  (Figure  29b),  there  is  a  small  scale  partial  mixing  that  occurs 
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adjacent  to  the  screen  that  serves  to  partially  mix  the  inflowing  groundwater  with  the  pre¬ 
pumping  well  water. 

As  mentioned  earlier,  the  transparent  face  plate  presents  a  no  slip  boundary  to  flow  within  the 
Half  Well.  The  Full  3D  well  avoids  the  artifacts  of  that  no-slip  boundary  and  allows  flow  to  be 
viewed  from  above  or  below  (where  the  screen  is  sealed  to  the  transparent  bottom  of  the  tank).  It 
is  more  difficult,  however,  to  produce  detailed  images  of  flow  within  the  Full  3D  Well.  Instead, 
arrays  of  thermistors  (i.e.,  monitoring  temperature)  and  E.C.  probes  (i.e.,  monitoring 
concentrations  of  NaCl  tracer)  were  installed  within  the  well.  These  probes  were  used  to  plot  the 
evolution  of  both  temperature  and  concentration  distributions  within  the  well. 

The  thermistor  string 
consisted  of  1.5mm  bead 
type  thermistors  with  a 
response  time  in  the  range 
of  a  few  seconds  in  order 
to  monitor  rapid  changes 
in  water  temperature  and 
minimize  disturbance  to 
flow.  These  thermistors 
were  calibrated  to  track 
each  other  within  0.005°C 
in  order  to  detect  the 
onset  of  unstable  thermal 
gradients.  Figure  35 
shows  the  temperatures 
within  the  Half  Well  as  a  function  of  depth.  Note  that  the  left-side  y-axis  relates  to  the  wellbore 
temperatures  and  the  small-scale  rapid  fluctuations  are  on  the  order  of  0.005°C.  The  right-hand 
y-axis  relates  to  the  ambient  temperatures  of  the  transparent  face  plate  of  the  Half  Well.  At  first 
(before  the  50  min.  mark)  as  the  ambient  temperatures  increased,  the  top  of  the  well  was  warmer 
than  the  bottom,  the  well  water  temperatures  were  uniform  and  warmed  simultaneously.  The 
uniform  well-water  temperatures  indicate  a  stable  thermal  regime,  i.e.,  warm  less  dense  water 
over  cooler  more  dense  water  (AT/Az  =  -l-0.10°C/m). 

Once  the  ambient  temperatures  reverse,  after  51  minutes,  the  well  water  thermal  gradients 
become  unstable,  i.e.,  cooling  upward  (AT/Az  =  -0.23°C/m).  Thermal  convection  is  indicated  by 
the  divergence  of  the  temperatures  and  the  fluctuation  of  temperature  with  time.  The  thermistor 
near  the  bottom  of  the  well  shows  a  periodic  fluctuation  observe  during  thermal  convection. 
These  indicate  small  plumes  of  cool  water  sinking  past  the  probe  from  farther  up  the  well.  The 
magnitude  and  period  of  these  fluctuations  are  used  to  estimate  the  rate  of  descent,  e.g., 
fluctuations  between  55  and  67  seconds  are  on  the  order  of  0.01 1°C  every  1.3  minutes  or  0.0085 


40 


50 


60 

Time  (min) 


70 


20 


19 


80 


U 

O 

5 

fO 

a 

E 

.V 


18 


2  Ambient  Temp. 


17 


Well  water 
temperature  at 
various  depths 


-110cm 

-80cm 

-30cm 

-20cm 

10cm 


- Well  Top 

—  Well  Bottom 


Figure  35.  Transition  from  thermal  stability  to  thermal  instability  with  the 
physical  model  due  to  inversion  of  ambient  temperature  gradients. 
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°C/min.  Dividing  that  temperature  fluctuation  by  the  thermal  gradient  (AT/Az  =  -0.23°C/m) 
gives  the  distance  that  water  traveled  and  the  velocity  of  descent,  (0.0085  “C/min)/(- 
0.23°C/m)=0.037m/min=3.7cm/min. 

When  the  onset  of  unstable 
thermal  gradients  is  more 
gradual  the  transition  to 
wellbore  convection  will  be 
more  gradual,  however,  the 
temperature  differentials  and 
fluctuations  still  occur  and  can 
be  used  to  characterize  the 
thermal  convection.  This  more 
gradual  transition  is  shown  in 
Figure  36  for  a  thermistor  20cm 
from  the  bottom  of  the  Full  3D 
Well.  The  vertical  gradients 
were  not  recorded  here  but  the 
transition  from  gradual 
convection  (before  120  min)  to  more  rapid  convection  is  evident.  This  indicates  that  the  rates 
and  vigor  of  thermal  convection  can  be  characterized  in  the  field  with  inexpensive,  commonly 
available,  single  thermistor  probes. 

The  electrical  conductance  (EC)  and  relative  concentration  (Cr)  of  NaCl  tracer  within  the  well 
models  is  measured  using  EC  probes  that  consist  of  two  parallel  8mm  long  stainless  steel  wires 
spaced  1mm  apart.  By  applying  a  IkHz  alternating  current  (1),  to  avoid  polarization,  and 
measuring  the  voltage  (V)  the  electrical  conductance  is  calculated  as  EC=EV.  With  a  reference 
solution  of  known  concentration  (Co)  and  measured  electrical  conductance  (ECq),  the  relative 
concentration  is  calculated  as  Cr=EC/ECo. 

The  initial  conditions  of  the  test  depicted  by  Eigure  29a  were  reproduced  in  with  NaCl  tracer. 
Before  pumping,  the  well  was  spiked  with  NaCl  tracer  and  mixed  to  give  a  uniform  Cr  =  1.0 
throughout  the  well.  Cr  was  then  monitored  at  various  levels  during  pumping  (at  1 .2  L/min)  to 
give  the  change  in  concentration  with  volume  pumped  (Nwv  on  Eigure  37)  The  concentration 
farthest  from  the  pump  intake  (the  top  of  the  well)  showed  the  slowest  decrease  in  concentration 
corresponding  to  the  zone  of  slow  moving  groundwater  in  the  analytical  and  numerical  models 
(Eigure  29b).  The  probe  at  the  bottom  of  the  well  (Eigure  37,  light  green  curve)  most  closely 
matched  the  concentrations  entering  the  pump  intake  (Eigure  37,  dark  blue  curve).  The  pumped 
concentrations  are  somewhat  lower  than  the  dashed  curve  which  indicates  the  theoretical 
concentrations  assuming  inviscid  flow  (Figure  30,  curve  A).  The  effect  of  dispersive  mixing 
(due  to  small  scale  eddy,  thermal  and  solutal  mixing)  is  to  cause  the  pre -pumping  well  water  to 
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Figure  36.  Temperature  fluctuations  indicating  thermal  convection  and 
thermal  instability  in  the  Full  3D  Well. 
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be  removed  more  rapidly  as 
shown  by  the  numerical  model 
of  a  pumped  well  (Figure  30, 
curve  E).  The  rapidly 
decreasing  concentration  at  the 
pump  (Figure  37,  dark  blue 
curve)  suggests  that  some 
degree  of  dispersive  mixing  is 
influencing  the  removal  of 
well  water  at  least  early  in 
pumping  (Nwv<2). 

The  physical  modeling 
suggests  that  natural  gradient 
flow  through  (horizontal  and 
vertical)  is  on  the  order  of 
millimeters  to  centimeters  per  day,  while  thermal  convection  during  unstable  conditions  causes 
movement  at  rates  on  the  order  of  centimeters  to  meters  per  hour  and  velocities  due  to  low-flow 
pumping  are  on  the  order  of  centimeters  per  minute.  Thus,  in-well  transport  via  free  convection 
is  generally  significant  only  in  the  absence  of  pumping  or  other  mechanisms  that  force 
convection.  Our  numerical  simulations  examining  in-well  transport  behavior  therefore  included 
thermal  and  solutal  convection  effects  in  studies  of  natural  gradient  groundwater  flow,  but 
neglected  those  effects  in  examination  of  in-well  transport  under  pumping  scenarios. 

5.3.2  Fracture  Borehole  Model 

The  physical  model  of  a  borehole  within  fractured  rock  was  designed  to  investigate  the 
interactions  between  fracture  flow  and  borehole  flow  under  controlled  conditions.  Imaging  the 
dye  passing  from  the  fracture  into  the  borehole  allowed  visualization  of  the  mixing  as  the  water 
entered  the  borehole  and  observation  of  mixing  within  the  borehole  (Figure  38).  The  zone  down 
gradient  of  the  borehole  represented  the  shadow  zone  indicating  the  influence  of  borehole  mixing 
on  the  groundwater  in  the  fracture.  The  lower  constant-head  reservoir  of  the  angled  fracture  was 
held  at  a  head  that  was  0.24cm  higher  than  the  head  of  the  upper  constant  head  reservoir.  This 
head  drop  resulted  in  a  hydraulic  gradient  of  0.00814,  a  volumetric  flow  rate  of  6.00cm  /min  and 
an  average  linear  velocity  of  5.00cm/min  toward  the  upper  constant  head  reservoir.  The  flow 
field  around  the  borehole  was  traced  with  Rhodamine  dye  introduced  as  a  uniform  constant 
source  in  the  high  head  reservoir  as  shown  at  the  bottom  of  Figure  38(a).  The  flow  field 
delineated  in  Figure  38(a)  shows  convergent  flow  up-flow  of  the  borehole  and  divergent  flow 
down-flow  of  the  borehole  as  shown  in  the  top  of  Figure  38(a). 


Number  of  Well  Volumes  Pumped  (N^,„) 


Figure  37.  Relative  concentration  changes  due  to  pumping  at  the  well 
bottom. 


The  dye  was  carried  into  the 
borehole  where  it  mixed  with 
the  borehole  water  and  was 
carried  out  the  down-flow 
side  of  the  borehole.  In  this 
case  the  dye  was 
approximately  0.01%  more 
dense  than  the  borehole 
water  and  mixed  by 
convection  (Figure  38(b) 
below  angled  fracture).  At 
room  temperature  a  0.01% 
increase  in  density  would  be 
equivalent  to  a  0.5°C 
decrease  in  temperature. 

Since  the  temperature  was 
not  controlled  in  this 
experiment  some  of  the  dyed 
water  in  the  fracture  was 
warmed  enough  to  reduce  the  density  below  that  of  the  borehole  water  and  rose  to  the  top  portion 
of  the  borehole  (Figure  38(b)  above  angled  fracture).  Because  the  dyed  water  mixed  with  the 
undyed  borehole  water  and  was  diluted,  a  zone  of  low  concentration  extended  down-flow  of  the 
borehole  producing  a  “shadow  zone”  of  borehole  water  as  shown  by  the  dilute  water  at  the  top  of 
Figure  38(a). 

For  the  next  series  of  experiments  the  concentration 
of  dye  was  reduced  by  a  factor  of  10  giving  a 
density  contrast  of  about  0.001%.  This  reduced 
density  contrast  reduced  the  sinking  velocity  and 
resulted  in  some  of  the  fracture  water  being  carried 
directly  across  the  borehole  before  mixing  with  the 
borehole  water.  Because  the  sinking  velocity  in  this 
experiment  was  on  the  order  of  millimeters  per 
minute  and  the  time  water  takes  to  traverse  was  less 
than  a  minute,  some  of  the  fracture  water  was 
carried  across  the  left  14  of  the  borehole  as  shown  in 
Figure  39.  This  indicates  that  the  degree  and 
distribution  of  mixing  between  fracture  water  and 
borehole  water  depends  on  the  balance  between 
fracture  flow  velocity  and  density  contrast.  In  this 
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High  Head,  Dye  Source 
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Figure  38:  Fracture  borehole  model  flow-field/borehole  interactions,  (a) 
downward  oblique  view  of  the  top  angled  fracture,  (b)  a  side  view  of  the 
fracture  model  showing  the  vertical  borehole,  the  angled  top. 
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case  the  fracture  water  entering  the  side  of  the  borehole  was  carried  across  and  exited  to  the 
fracture  while  the  water  in  the  center  of  the  borehole  partially  mixed  with  the  borehole  water. 
This  series  of  experiments  indicate  that  fracture  flow  into  the  borehole  and  mixing  within  the 
borehole  may  result  in  heterogeneous  concentration  distributions  along  the  borehole  and  a 
shadow  zone  of  altered  water  up-gradient  of  the  borehole. 

5.5.5  Planar-laser-induced  fluorescence  (PLIF)  tomographic  imaging 

The  PLIF  tomography  system  developed  successfully  demonstrated  the  ability  to  record  dye 
movement  in  a  pseudo  2D  model  well  into  3D  data  describing  concentration  changes  over  time 
(Figure  40).  The  2D  and  3D  imagery  developed  from  that  effort  clearly  illustrated  the  3D  nature 
of  the  dye  transport  within  the  well,  including  the  initial  entry  point,  and  the  increasing 
complexity  of  the  dye  plumes  over  time.  However,  while  the  imaging  effort  allowed 
development  of  a  novel  interrogation  method,  and  a  dataset  that  helped  to  illustrate  the 
complexity  of  in-well  flow  behavior,  we  concluded  that  simpler  non-physical  analysis  methods 
provided  more  efficient  means  of  testing  hypotheses  about  in-well  mixing. 


Figure  40.  Single-frame  of  PLIF  imagery  (A)  and  3-D  iso-surface  reconstructions  (B)  based  on  ~30  frames  of  video 
during  laser  translation  through  well.  Horizontal  flux  during  experiment  was  =  160  cm  day  ’.;  injection  consisted  of 
approximately  20  ml  of  100  ppm  rhodamine  WT  dye.  Iso-surfaces  represent  a  constant,  arbitrary,  concentration 
level  measured  via  fluorescence,  not  calibrated  to  actual  dye  concentration. 
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5,4  Physical  and  Numerical  Modeling 

5.4. 1  Non-Thermal  Density  Driven  Flow 

A  Rayleigh  number  analysis  strongly  suggests  that  natural  conveetion  should  oecur  in  the  open 
water  section  of  Britt’s  (2005)  apparatus  given  sufficient  density  (including  temperature) 
contrast.  To  examine  free  convection  in  that  system,  in  which  horizontal  flow  should  act  to 
attenuate  the  vertical  extent  of  natural  convective  transport,  we  simulated  flow  and  transport  in 
the  domain  depicted 
schematically  in  Figure  41. 

The  2D  domain  essentially 
reproduces  Britt’s  experiment 
in  dimensions, 

hydrostratigraphy,  injection 
behavior  and  fluid  flow  rates. 

Fluid  flow  was  controlled  at 
the  left  and  right  boundaries 
of  the  domain,  with  a 
specified  horizontal  flux  of 
6E-7  m/s  (5.2  cm/day)  at  the 
left-hand  domain  boundary 
and  an  arbitrary  specified 
head  at  the  right-hand 
boundary.  The  porosity  of  the 
porous  media  is  universally 
defined  as  30%,  so  that  the 
seepage  velocity  in  the  Darcy’s  law  domains  of  the  model  closely  matches  that  described  by  Britt 
(2005).  Boundary  conditions  for  solute  transport  included  a  specified  zero-concentration 
boundary  on  the  left-hand  side  and  an  advective  flux  boundary  on  the  right.  A  source  point  term 
was  defined  1 5  cm  upstream  of  the  well,  centered  vertically  within  the  domain,  with  tracer  as  in 
Britt’s  experiment  -  a  solution  injection  rate  equal  to  1/3  the  total  inflow  rate  of  the  upstream 
boundary.  A  permeability  of  9.1E-11  m  was  specified  for  the  aquifer  and  the  permeability  of 
the  sand  pack  adjacent  to  the  well  was  9.1E-10  m  .  With  the  configuration  described, 
isopotentials  diverge  near  the  injection  point  but  are  essentially  parallel  when  they  reach  the  well 
(Eigure  42).  Only  minor  deflection  of  the  flow  lines  occurs  in  the  well  section  of  the  model. 
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In  the  initial  simulation,  the 
density  of  the  injeeted 
solute  was  set  to  match  that 
of  the  background  solution 
and  in  that  case  the  density 
behaves  as  one  might 
intuitively  expect,  being 
transported  across  the  well 
with  only  minor  dispersive 
behavior  due  primarily  to 
molecular  diffusion  (Figure 
43a).  The  effect  of  solute 
plume  density  on  transport 
through  the  well  was  then 
examined  by  repeatedly 
running  the  simulation  with 
increasingly  higher  injection 
density.  Natural  convection 
began  at  an  approximate  relative  density  increase  of  lE-8.  At  that  point,  the  competition 
between  natural  convection  and  forced  convection  creates  extensive  mixing  in  the  well  (Figure 
43b),  but  only  in  the  downward  direction,  consistent  with  the  net  buoyant  force  on  the  solution. 
The  density  driven  flow  evidences  the  characteristic  finger-like  projections  of  Rayleigh  Taylor 
instabilities  (Bird  et  al.  1960),  with  mushroom-like  caps  at  the  distal  end  of  the  downward 
directed  plume.  Lower  density  water  entering  below  the  solute  plume  also  forms  finger-like 
projections  that  rise  upward  through  the  solute  plume,  until  reaching  the  level  of  clean  water 
entering  above  the  solute  plume.  The  competition  between  the  descending  plume  and  lower 
density  water  rising  through  it  forces  the  projections  of  each  movement  to  change  position  with 
time,  forming  a  sort  of  wagging  tail  in  which  the  translation  is  proportional  to  the  density 
contrast  of  the  entering  fluids. 
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Figure  42.  Isopleths  of  total  water  potential  describing  the  flow  regime  in 
simulations  of  Britt’s  2005  sand  tank  experiments.  Streamlines  (red  curves) 
and  flow  vectors  illustrate  minor  deflection  in  the  well  section  of  the 
domain. 
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Injection  density  =  background  density  I  Injection  density  = +lE-8po 


Figure  43.  Simulated  solute  transport  in  a  numerical  replication  of  Britt’s  2005  experiment  with  (A)  no  increase  in 
solution  density  and  (B)  a  +1E-8  increase  in  solution  density  relative  to  background  density.  Panel  C  illustrates  the 
effect  of  adding,  to  the  simulation  shown  in  panel  B,  a  partial  barrier  (80%  of  well  width)  in  the  well  section  of  the 
domain.  Velocity  vectors  (black  arrows)  are  shown  in  well  and  sand -pack  domains.  Contours  are  total  water 
potential. 


As  indicated  previously,  the  slight  inerease  in  density  producing  the  observed  solutal  eonvection 
reflects  a  solute  density  eontrast  of  only  ~10  |ug/L.  Consistent  with  the  Rayleigh  analysis,  the 
simulation  emphasizes  that  very  minor  solute  concentrations  ean  exert  profound  effects  on 
vertical  distribution  of  solute  in  a  well.  This  implies  that  in  situ  monitoring  devices  would  in 
many  oases  measure  a  looal  ooncentration  that  not  refleot  either  the  adjacent  formation 
concentration  or  a  flow-weighted  average  of  well  inflow  oonoentrations.  This  is  generally 
undesirable  behavior  if  one  of  the  goals  of  using  in  situ  monitoring  devices  is  to  measure  the 
vertical  distribution  of  the  formation  or  an  average  of  the  whole  screen  seotion.  ProHydro,  Ino., 
manufacturer  of  the  Snap  Sampler — an  in  situ  sampling  devioe  that  is  activated  remotely — has 
experimented  with  use  of  in-well  flow  inhibitors  to  attempt  to  isolate  separate  sampling  deviees 
in  an  open  well  sereen  section  (Figure  44).  The  deviees  do  not  give  eomplete  hydrologic 
isolation  but  can  provide  inereased  stratifieation  of  eoneentration  measurements  from  multi-level 
sampling.  Because  some  of  the  vertieal  movement  induced  by  natural  conveetion  reflects 
conservation  of  momentum,  it  is  possible  that  even  an  incomplete  barrier  within  the  well  could 
significantly  inhibit  vertical  convection.  As  a  test  of  this  hypothesis,  we  plaeed  a  partial  barrier 
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Figure  44.  In  well  flow  inhibitor  used  by 
ProHydro,  Inc.  to  limit  vertical  solute 
migration  in  a  well. 


in  the  well  section  domain  of  our  numerical  model  of 
Britt’s  (2005)  apparatus.  The  width  of  the  barrier  is  80% 
of  the  width  of  the  well  section.  With  this  partial  barrier  in 
place,  the  downward  convective  flow  induced  by  the 
slightly  higher  density  of  the  injected  solution  is  largely 
eliminated  at  the  location  of  the  barrier.  While  this 
simulation  represents  a  2D  flow  regime,  and  cannot  be 
considered  to  accurately  represent  3D  flow  in  a  well,  the 
simulation  suggests  that  partial  barriers  may  provide  an 
effective  means  of  isolating  in  situ  sampling  devices  under 
certain  conditions. 

The  simulations  thus  far  presented  demonstrate  that  very 
slight  concentration  differences  between  the  injected 
solution  and  the  background  solution  can  induce 
significant  vertical  convection.  In  Britt’s  experiment 
however,  the  injected  dye  was  observed  to  migrate  rapidly 
both  upward  and  downward  in  the  well,  whether  the 
injection  solution  density  was  slightly  higher  or  lower  than  background  density.  In  the  course  of 
investigating  that  phenomenon,  INL  constructed  a  similar  sand  tank  and  conducted  several 
preliminary  experiments  using  a  rhodamine  WT  solution  with  density  approximately  1.0001  x 
that  of  the  background  solution.  In  each  experiment  conducted  with  the  tank  set  up  similar  to  the 
arrangement  of  Britt’s 
2005  experiment,  dye 
entering  the  well  was 
observed  to  migrate 
faster  in  the  vertical  than 
in  the  horizontal 
direction,  even  at 
horizontal  velocities  up 
to  2.5  m/d  and  with  a 
higher  permeability  layer 
arranged  in  the  middle  of 
the  tank  to  attempt  to 
reduce  vertical  transport 
within  aquifer  portions 
of  the  tank.  A  useful 
effect  of  the  higher 
velocity  experiments 
was  greater  separation  of 
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Figure  45.  Results  of  a  dye  tracer  experiment  conducted  at  INL  investigating 
transport  to  a  well  of  a  negatively  buoyant  injection  solution.  The  upward 
movement  of  the  dye  upon  entry  into  the  well  is  hypothesized  to  reflect  thermal 
convection  due  inadvertent  warming  of  the  injection  solution. 
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horizontal  and  vertical  movement  within  the  well.  In  the  experiment  illustrated  in  Figure  45,  for 
example,  the  injected  solution  clearly  rose  rapidly  upon  entering  the  well  and  then  descended 
after  traversing  a  distance  of  approximately  half  the  well  diameter  (Figure  46). 

Subsequent  discussions  about  this  phenomenon  indicated  that  nearly  identical  behavior  had  been 
observed  in  a  similar  experiment  at  the  University  of  Toledo.  Because  in  each  case,  the  injected 
solutions  were  believed  to  be  negatively  buoyant  based  on  solute  concentration  and  density,  we 
hypothesized  that  inadvertent  temperature  differential  (warming)  of  the  injected  solution  might 

be  causing  thermal  convection  in 
the  well. 
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Figure  46.  Results  of  a  dye  tracer  experiment  conducted  at  INL 
investigating  transport  from  an  aquifer  to  a  well  with  buoyant  effects 
operating  in  the  short  dimension  of  the  tank  (out  of  the  page  in 
figure).  In  this  experiment,  vertical  dye  transport  was  notably  more 
diffuse  and  slower  than  in  comparable  experiments  with  the  tank 
arranged  in  vertical  orientation. 


To  test  the  hypothesis  that  rapid 
vertical  migration  of  the  dye  within 
the  well  sections  of  the  tank 
resulted  from  buoyant  effects,  we 
conducted  a  single  similar 
experiment  with  the  long  axis  of  the 
well  in  the  horizontal  plane.  With  a 
horizontal  water  flux  of 
approximately  7E-6  m/s  (2  ft/day), 
the  dye  entering  the  tank  behaved 
much  more  as  would  be  expected 
for  diffusive  transport  in  a  laminar 
flow  regime.  Unfortunately, 
increasing  water  leakage,  and  air 
entry,  from  the  upper  seals  in  the 
assembly  limited  our  ability  to 
continue  such  experiments. 
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To  examine  the  potential  for  the  combined  effects  of  thermal  and  solutal  convection  in  these 
sand  tank  experiments  we  added  heat 
transport  to  our  numerical  model  of 
Britt’s  2D  sand  tank  experiment. 

Results  of  a  simulation  in  which  the 
injected  solution  was  given  a  constant 
injection  temperature  of  +0.001°C 
and  a  relative  solutal  density  increase 
of  lE-8  are  shown  in  Figure  47. 

Simulated  solution  behavior  in  the 
well  section  of  the  model  is  observed 
to  have  many  of  the  characteristics 
displayed  in  our  preliminary 
laboratory  results  (Figure  45).  Flpon 
entering  the  well,  the  dye 
immediately  rises  due  to  thermal 
convection  and  subsequently 
descends  as  the  injected  solution 
mixes  with  and  conducts  heat  to  the  surrounding  cooler  water  and  in  both  simulated  and 
observed  results  the  height  of  that  rise  decreases  over  time.  The  simulation  also  appears  to 
effectively  reproduce  the  behavior  of  the  injection  solution  below  the  entering  plume,  where  the 
descending  plume  displays  the  characteristic  lingering  of  Rayleigh-Taylor  instabilities.  While 
the  comparison  described  here  is  preliminary  and  qualitative,  the  magnitude  of  the  thermo- 
solutal  convection  effects  are  consistent  with  Rayleigh  number  analysis.  We  take  these  results  as 
strong  evidence  that  our  simulations  effectively  capture  the  combined  effects  of  thermal  and 
solutal  convection  in  these  2D  representations  of  a  well.  Subsequent  simulation  studies  will 
attempt  to  examine  these  effects  in  a  3D  simulation  of  a  well. 

The  vertical  extent  of  density-induced  mixing  clearly  depends  on  the  density  contrast,  the  type  of 
convection  present  (thermal,  solutal,  or  both)  and  the  horizontal  flux  across  the  well.  The 
interaction  between  horizontal  flux  and  vertical  convection  is  difficult  to  assess  without  detailed 
analysis  of  the  flow  fields.  As  an  example.  Figure  48  illustrates  how  a  doubling  of  the  injection 
density  affects  vertical  transport  in  a  2D  model  similar  to  that  used  to  reproduce  Britt’s  (2005) 
experiment  but  extended  in  the  vertical  dimension  to  examine  potential  differences  in  the  extent 
of  vertical  mixing.  With  a  5E-9  relative  increase  in  density,  the  convective  mixing  effect  appears 
to  extend  approximately  30  cm  below  where  the  plume  enters  the  well.  When  the  relative 
increase  is  doubled,  to  lE-8,  the  vertical  extent  of  mixing  is  not  doubled,  but  extends 
approximately  another  20  cm  below  the  former  case.  Doubling  the  relative  density  again  (Figure 
48C),  however,  causes  the  vertical  mixing  effect  to  extend  to  the  bottom  of  the  well, 
approximately  30  cm  further  than  in  Figure  48B. 


Figure  47.  Results  of  a  numerical  simulation  of  Britt’s  (2005) 
experiment  with  a  concentration-based  relative  density  increase  of 
lE-8  and  a  temperature  0.00 1°C  above  background.  Velocity 
vectors  (black  arrows)  are  shown  in  well  and  sand-pack  domains. 
Contours  are  total  water  potential. 
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In  the  essentially  2D  flow  regime  of  Britt’s  (2005)  sand  tank  configuration,  the  effects  of  the 
descending  solute  plume  extend  to  approximately  a  meter  with  a  relative  density  difference  of 
5E-8  and  a  horizontal  flow  velocity  of  5.2  cm  day'\ 


Injection  density  =  +5E-9po  Injection  density  =  +lE-8po  Injection  density  =  +2E-8po 

q  =  5.184  cm/day  q  =  5.184  cm/day  q  =  5.184  cm/day 
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Figure  48.  Simulations  designed  to  examine  how  the  vertical  extent  of  convective  mixing  varies  with  density 
contrast.  In  both  simulations  the  horizontal  flux  is  6E-7  m/s  (5.2  cm/day).  In  panel  A,  the  injection  density  is  5E-9 
above  that  of  the  background  solution.  In  panels  B  and  C,  the  relative  density  increase  is  double  (B)  and  quadruple 
(C)  that  of  A. 

5.4.2  Thermal  convection  in  a  well 

Although  numerous  studies  have  examined  the  critical  Rayleigh  number  for  thermal  convection 
in  a  well  (eg.  Hales,  1937;  Vanaparthy  et  ah,  2003),  the  stability  analysis  does  not  characterize 
the  convective  movement,  and  a  description  of  the  magnitude  and  scale  of  the  convective 
movement  is  needed  to  understand  how  it  might  affect  concentration  distribution  in  a  well. 
Other  studies  have  examined  the  convective  processes  by  measuring  in-well  temperatures 
(Gretener,  1967;  Diment,  1967;  Cermak  et  ah,  2007),  but  those  data  provide  information  about  a 
limited  number  of  scenarios,  and  likely  include  other  effects  that  complicate  interpretation  of  the 
data.  The  thermal  gradients  that  drive  convection  close  to  the  critical  Rayleigh  number  are  also 
sufficiently  small  that  the  temperature  variations  that  describe  the  convective  movement  are 
difficult  to  measure  accurately.  Recently,  Berthold  and  Bomer  (2008)  conducted  transient  2D 
simulations  of  thermally  convective  Navier  Stokes  flow  to  compare  with  observations  of  laser 
sheet  illuminated  convection  in  a  tube  in  the  laboratory.  Like  the  above  simulations  of  mixed 
forced  and  free  convection  however,  those  simulations  describe  a  2D  rectangular  flow  domain. 
To  examine  convective  behavior  under  conditions  typical  in  shallow  wells  we  conducted  2D 
axisymmetric  flow  simulations  designed  to  reproduce  observations  made  by  Martin-Hayden  and 
Britt  (2006)  of  thermal  convection  in  instrumented  unscreened  well. 
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Martin-Hay  den  installed  three  strings  of  high-resolution  thermistors  in  and  adjacent  to  an 
initially  empty  2”-diameter  PVC  well  with  a  bottom  depth  of  5.7  meters.  One  thermistor  string 
monitored  the  vertical  temperature  profile  just  outside  the  well,  another  monitored  temperature 
along  the  inside  wall  of  the  well  and  another  monitored  temperature  along  the  axis  of  the  well. 
After  filling  the  tube  with  water  to  a  depth  of  2.2  meters,  temperatures  were  monitored  to 
observe  behavior  while  and  following  equilibration  with  the  background  temperature  gradient. 
Our  numerical  simulation  reproduced  that  experiment  by  incorporating  a  non-flowing  (thermal 
transport  by  diffusion,  no  advection)  domain  around  a  2D  axisymmetric  well  section  (Figure 
49A)  in  which  Navier-Stokes  equations  are  solved.  Dimensions,  temperatures  applied  at  top  and 
bottom  boundaries,  and  temperature  gradient  applied  at  the  right-hand  boundary  of  the  model 
(~l°C/m)  were  defined  to  reproduce  the  field  experiment  of  Martin-Hayden  (Martin-Hayden  and 
Britt,  2006).  Results  of  this  simulation  demonstrate  a  toroidal  convection  pattern  in  the  well  with 
two  separate  cells,  with  upward  convection  (Figure  49B)  deformation  of  that  gradient  close  to 
the  well  because  of  the  toroidal  convection  pattern  in  the  well.  Colors  and  arrows  in  the  well 
show  intensity  and  direction  of  convection  in  the  well,  with  a  maximum  velocity  of  -8  cm/min. 
The  results  reproduce  many  of  the  features  of  Martin-Hayden ’s  experiment,  including  the 
apparent  number  of  convection  cells,  estimated  vertical  velocity,  and  the  effect  of  convection  on 


Figure  49.  Schematic  diagram  (A)  and  results  (B)  of  a  2D  axisymmetric  model  of  thermally  convective  flow  in  a 
well  in  the  shallow  subsurface.  Dimensions,  temperatures  applied  at  top  and  bottom  boundaries,  and  temperature 
gradient  applied  at  right-hand  boundary  were  designed  to  reproduce  the  field  experiment  of  Martin-Hayden  and  Britt 
(2006).  Isotherms  (B)  show  background  gradient  of  approximately  l°C/m  and  deformation  of  that  gradient  close  to 
the  well  because  of  the  toroidal  convection  pattern  in  the  well.  Colors  and  arrows  in  well  show  intensity  and 
direction  of  convection  in  the  well. 
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vertical  temperature  profiles  (Figure  50)  in  the  well.  Simulated  temperature  profiles  domain 
demonstrate  crossover  between  the  inner  and  outer  profiles  that  is  similar,  though  greater,  to  the 
crossover  observed  between  the  inner  and  outer  thermistor  strings.  Differences  between 
simulated  and  observed  temperature  distribution  may  be  due  the  representation  of  thermal 
conductivity  immediately  adjacent  to  the  well,  as  the  simulated  temperatures  adjacent  to  the  well 
clearly  demonstrate  stronger  connection  to  in-well  temperatures  than  Martin-Hayden  observed. 

To  our  knowledge,  this  is  the  first  study  of  thermal  convection  in  a  shallow  well  that  combines 
numerical  simulation  with  precise  temperature  monitoring  at  several  radial  positions  in  and 
around  the  well,  and  for  which  simulation  results  appear  to  provide  a  reasonable  characterization 
of  the  observed  convective  effects.  Additional  simulation  experiments  are  currently  being 
performed  to  examine  several  related  processes,  including  (1)  the  rate  of  vertical  mixing  induced 
by  thermal  gradients  with  and  without  solutal  convective  effects,  (2)  the  likely  mechanisms  that 
produce  in-well  temperature  oscillations  that  are  frequently  observed  in  in-well  convection 
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Figure  50.  Temperatures  measured  in  a  shallow  unscreened  well  (Martin-Hayden  2001)  and  calculated 
temperatures  (B)  in  a  numerical  simulation  designed  to  reproduce  that  experiment. 
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studies  and  (3)  convective  effects  in  multilayered  formations  with  varying  concentration, 
permeability  and  total  water  potential. 

Numerical  simulations  have  thus  far  been  used  to  investigate  the  mechanisms  causing  the 
counter  intuitive  in-well  solute  behavior  described  by  Britt  (2005),  examine  the  character  of 
solutal  and  thermal  convection  with  and  without  competing  forced  convection  perpendicular  to 
the  free  convection  gradient,  and  to  examine  in-well  solute  concentration  during  low-flow 
pumping.  In  each  of  these  studies,  good  matches  between  observed  and  simulated  behavior 
appears  to  confirm  that  solution  of  the  Navier- Stokes  equations  for  in-well  flow  is  critical  for 
accurate  representation  and  prediction  of  in-well  mixing  effects  under  ambient  flow  and,  in  some 
cases,  low-flow  pumping  conditions.  Results  further  demonstrate  that  solutal  and  thermal 
convective  effects  can  produce  substantial  vertical  mixing  under  ambient  flow  conditions  in  a 
well  and  that  these  effects  must  be  considered  in  well  design  where  non-purge 
sampling/monitoring  methods  are  employed.  Finally,  simulations  of  low- flow  pumping  in  a  well 
are  consistent  with  analytical  solutions  describing  temporal  changes  in  effluent  concentration 
during  pumping.  Those  calculations  demonstrate  that  (1)  the  volume  that  must  be  purged  from  a 
well  to  ensure  a  flow-weighted  average  of  the  surrounding  formation  concentrations  should  be 
based  on  hydraulic  considerations  rather  than  purge  parameters  that  may  suggest  stability  before 
the  proper  average  concentration  is  reached  and  (2)  while  the  necessary  purged  volume  may  vary 
with  the  degree  of  in-well  mixing  induced  by  irregularities  or  other  perturbations  from  the  ideal 
well  systems  considered  here,  under  ideal  conditions  removal  of  three  well  volumes  is  sufficient 
to  reach  the  desired  flow-weighted  average  of  the  concentration  distribution  outside  the  well. 
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6.0  Conclusions  and  Implications 

Physical  and  numerical  modeling  efforts  were  completed  to  examine  the  relationship  between 
solute  distribution  in  a  well  and  that  in  the  adjaeent  formation.  Physieal  mixing  processes  that 
oeeur  during  ambient  flow  in  aquifers  and  during  aetive  pumping  aetivities  were  explored  to 
show  how  potentially  stratified  eontaminants  may  manifest  in  open  well  sereens.  Complex  flow 
patterns  driven  by  flow  velocity  changes,  very  small  density  contrasts,  and  temperature  eontrasts 
were  identified  as  possible  sourees  of  in-well  eonvection  and  mixing. 

These  varied  drivers  for  eontaminant  redistribution  work  together  to  alter  in-well  coneentration 
distributions  from  that  in  the  adjaeent  formation.  These  effeets  may  lead  to  virtually  eomplete 
mixing  within  the  well  sereen.  This  is  not  to  say  that  all  wells  are  homogenized  with  respect  to 
aquifer  stratifieation,  but  they  do  refleet  some  degree  of  redistribution  that  is  nearly  always 
present.  Thus,  while  stratification  can  be  found  in  wells,  it  is  unlikely  that  sueh  stratifieation  is 
equivalent  to  that  found  in  the  adjacent  aquifer.  When  little  or  no  stratification  is  found  in  a 
monitoring  well,  it  is  likely  the  eoneentration  identified  throughout  the  well  represents  an  aquifer 
with  little  stratification  in  the  target  interval,  a  specific  interval  that  is  weighted  due  to  hydraulic 
effects,  or  an  averaging  of  the  screen  interval  exposed  to  the  aquifer.  Rather  than  a  problem,  this 
eould  be  viewed  as  an  advantageous  condition  where  purging  wells  probably  adds  no  significant 
benefit  for  data  utility.  This  investigation  did  not  attempt  to  quantify  how  frequent  that  condition 
exists,  but  others  have  found  that  it  is  quite  common  (ITRC  2007).  Therefore,  given  the  variety 
of  drivers  toward  strong  redistributive  effects,  the  default  assumption  should  be  that  a  given  well 
is  probably  mixed  unless  shown  otherwise.  In  the  case  where  pumped  and  passive  samples  do 
not  reasonably  mateh,  detailed  investigation  may  be  helpful  for  determining  the  cause.  However, 
determining  such  root  cause  must  be  weighed  against  cost,  and  may  not  be  worthwhile  unless 
eritieal  deeisions  are  being  made  on  a  well-speeifie  basis.  There  are  conditions  that  can  cause 
eoneentration  dependenee  on  vertical  position,  including  positive  vertieal  temperature  gradients, 
but  these  eonditions  are  transient  and  also  do  not  neeessarily  maintain  an  aquifer/well 
concentration  elevation  parity. 

The  major  eonclusions  developed  from  the  work  include  the  following; 

I.  Field  Setting 

A.  Identification  of  contaminant  stratification  in  the  aquifer  using  longer  sereen  wells 
(>5  ft)  is  dependent  on  multiple  faetors  including  contaminant  inflow  position, 
sampling  method,  and  degree  of  zone  isolation. 

1.  Contaminant  position  in  the  aquifer  relative  to  the  passive  sampling  deviee 
position  in  the  well  is  usually  not  known,  and  understanding  of  that  geometry 
may  not  be  substantially  improved  by  passive  sampling  in  an  open  bore. 
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2.  Similarly,  aquifer  contaminant  position  relative  to  the  pump  intake  position  is 
rarely  known,  but  may  be  discerned  with  multiple  samples  collected  during 
purging. 

3.  Isolated  zone  sampling  can  improve  determination  of  aquifer  contaminant 
stratification,  using  packers  or  other  mixing  inhibition  devices  installed  in- 
between  passive  samplers. 

B.  Representation  of  the  adjacent  aquifer  through  pump  sampling  is  dependent  on 
several  factors  including  contaminant  inflow  position,  volume  of  the  purge,  and 
contaminant  of  interest. 

1.  Contaminant  position  in  the  aquifer  adjacent  to  the  well  relative  to  the  pump 
position  drives  the  speed  to  which  contaminant  concentration  stability  is 
achieved. 

2.  Volume  of  water  removed  from  the  well  is  a  more  reliable  predictor  of 
contaminant  concentration  stability  than  measurement  of  traditional  “purge 
parameters”  such  as  temperature  pH,  EC,  ORP  and  DO;  however,  even  large 
volumes  removed  do  not  always  assure  contaminant  concentration  stability. 

3.  Some  contaminants  of  interest  stabilize  faster  than  others,  depending  on 
conditions  at  the  individual  well,  suggesting  that  contaminant  distribution 
heterogeneity  or  chemical-specific  biological  activity  can  influence  the 
stability  of  contaminant  concentrations  during  purging. 

C.  Sample  handling  at  surface  can  impart  substantial  artifacts  to  measured 
concentrations; 

1.  Pouring  samples  for  VOC  analysis  in  the  open  air  results  in  differential  losses 
dependent  on  the  properties  of  the  chemical  of  concern. 

2.  Conventional  thought  that  pouring  sample  “smoothly”  with  little  agitation  is 
an  important  aspect  to  VOC  recovery  is  largely  incorrect — rather,  surface  area 
of  exposure  is  the  apparent  driver  of  agitation-based  losses. 

3.  Conventional  thought  that  pouring  a  sample  slowly  is  an  important  aspect  of 
VOC  recovery  is  also  largely  incorrect — rapid  bottle  filling  with  short 
exposure  time  appears  to  limit  VOC  losses. 

4.  Pouring  and  sealing  sample  quickly  (at  250ml/min  or  greater — 10  seconds  or 
less  to  fill  a  vial)  and  filling  with  the  discharge  tube  submerged  (very  little 
surface  area  exposed)  yielded  up  to  a  combined  nearly  20%  higher  average 


ER-1704  Final  Report  2014 


VOC  recoveries  compared  to  pouring  slowly  (50ml/min)  and  down  the  side  of 
the  container. 

II.  Laboratory  Experimentation 

A.  Physical  model  experiments  confirm  previous  work  by  Britt  (2005)  that  small 
heterogeneities  exert  strong  influence  on  flow  in  the  open  bore  of  a  model  well. 

1.  Horizontal  flow  of  near  neutrally  buoyant  tracer  in  porous  media,  including 
both  simulated  aquifer  and  filter  pack,  is  confirmed,  given  a  horizontal  flow 
field. 

2.  Horizontal  flow  is  not  maintained  in  the  open  well  bore  upon  tracer  incidence 
into  the  non-porous  media,  for  all  densities  tested. 

3.  Similar  lack  of  horizontal  flow  and  in- well  mixing  phenomena  were  observed 
in  simulated  fractured  media. 

B.  Horizontal  laminar  flow  across  the  model  bore  could  not  be  produced  (water 
entering  one  side  of  the  model  well  exiting  at  the  same  vertical  interval  with 
limited  mixing). 

1.  No  configuration  of  velocity  or  lithologic  heterogeneity  could  reproduce 
horizontal  flow  in  the  open  bore. 

2.  Density  could  not  be  controlled  sufficiently  to  test  all  possible  scenarios — 
perfect  neutral  buoyancy  was  the  most  difficult  to  approach,  and  given  that 
difficulty,  we  conclude  that  it  is  equally  or  perhaps  even  more  unlikely  to 
occur  in  natural  environments. 

C.  Slight  vertical  hydraulic  gradients  in  a  model  aquifer  apparatus  yield  flow  in  the 
well  with  similar  trajectories  predicted  by  analytical  approaches;  namely,  very 
small  vertical  gradients  in  porous  media  are  amplified  by  the  high  hydraulic 
conductivity  of  the  well  bore  resulting  in  substantial  vertical  redistribution. 

1.  An  upward  vertical  gradient  of  0.005  m/m  (10%  of  the  0.05  m/m  horizontal 
gradient)  in  the  U  Toledo  tank  yielded  a  vertical  redistribution  of  dye  in  the 
well. 

2.  The  same  gradient  caused  dye  to  be  distributed  over  an  approximately  20-cm 
vertical  distance  down-gradient  of  the  model  well. 

D.  Slight  density  and/or  temperature  contrasts  are  effective  causes  of  tracer 
redistribution  in  the  aquifer/well  models. 
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1 .  All  tracer  density  contrasts  yielded  redistribution  in  the  well  models 

2.  In  examples  where  temperature  contrast  was  also  suspected,  compound 

behaviors  changing  density  of  the  tracer  were  observed. 

III.  Modeling  of  Field  and  Laboratory  Experiments 

A.  A  series  of  numerical  modeling  experiments  were  completed  in  which  solute 
transport  is  tied  to  solution  of  the  Navier-Stokes  equations  in  the  well  and  to 
Darcy’s  law  only  outside  the  well.  Good  matehes  between  simulation  results, 
experimental  data,  and  theoretical  analysis  of  flow  and  transport  support  our 
hypotheses  of  homogenization  effects. 

B.  Numerical  modeling  experiments  confirm  our  conceptual  model  of  in- well  flow 
and  transport  behavior  that  flow  is  best  represented  by  a  Poseuille  flow  model. 

i.  Simulations  of  well  effluent  concentrations  during  pumping  indicate  that 
purging  three  well  volumes  effeetively  brings  the  well  back  to  a  flow- 
weighted  average  of  eoneentrations  in  the  surrounding  formation. 

ii.  Simulations  of  well  effluent  concentrations  during  pumping,  using 
Darcy’s  law  to  describe  in-well  flow,  indieate  that  approximately  five  well 
volumes  would  be  required  to  reach  the  desired  flow-weighted  average. 

C.  Numerical  modeling  experiments  demonstrate  that  very  slight  density  differenees 
in  an  open  well,  due  to  solutal  or  thermal  effects,  can  induce  signifieant  vertieal 
convection  and,  therefore,  much  greater  vertical  mixing  than  would  otherwise  be 
anticipated. 

i.  Density  eontrasts  equivalent  to  as  little  as  10  parts  per  billion  of  dissolved 
solids  are  enough  to  eause  near  complete  vertical  redistribution  of  solute 
in  the  short  well  seetion  ineorporated  in  Britt’s  (2005)  apparatus. 

ii.  Zone  isolation  within  the  bore  using  an  80%  impermeable  barrier  was 
suffieient  to  effectively  stop  the  vertical  density-driven  flow. 

iii.  Convection  rates  during  mild  unstable  gradients,  e.g.,  l“C/m,  can  drive 
eonvection  with  velocities  on  the  order  of  10  cm/min. 

iv.  While  typical  geothermal  gradients  are  sufficient  to  produce  convection  in 
the  wells  of  diameter  typieal  of  observation  wells,  much  stronger 
temperature  gradients  near  the  ground  surface  where  the  seasonal 
temperature  signal  penetrates  to  a  depth  of  approximately  20  meters. 
Simulation  of  experimental  data  from  a  near  surface  well  in  this  study 
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demonstrates  that  we  can  reasonably  accurately  describe  convective 
effects  that  likely  exert  significant  control  on  measured  concentrations  in 
near-surface  monitoring  wells.  Further  efforts  will  attempt  to  characterize 
convective  behavior  and  its  effects  at  lower  and  higher  temperature 
gradients  than  described  in  this  report. 

V.  Examples  of  additional  effects  to  be  considered  include  the  differences  in 
convective  behavior  induced  by  seasonal  thermal  changes,  which 
redistribute  solutes  in  a  well  differently  at  different  times  of  the  year. 

1.  During  winter  the  cool  surface  temperatures  lead  to  unstable 
thermal  gradients  within  monitoring  wells,  i.e.,  cooler  denser  water 
over  warmer  less  dense  water  leads  to  thermal  convection. 

2.  As  surface  temperatures  warm,  stable  thermal  gradients  will 
propagate  downward  at  a  rate  determined  by  the  rate  of  warming, 
the  thermal  diffusivity  of  the  formation  and  groundwater  recharge. 

3.  Deeper  regions  where  geothermal  gradients  persist  (i.e.,  increasing 
temperature  with  depth)  will  be  continually  thermally  unstable. 


Overall  Conclusions  and  Technical  Transfer.  Contaminant  redistribution  effects  in  wells  are 
nearly  always  present.  Complete  mixing  appears  to  be  very  common;  however,  it  is  not 
universal.  There  is  a  continual  balance  between  inflowing  contaminant  stratification  (where 
present)  and  factors  driving  in-well  mixing.  Findings  here  imply  common  and  very  small  drivers 
are  responsible  for  slow  but  vigorous  mixing  relative  to  the  residence  time  of  water  flowing 
through  a  typical  well  screen.  Therefore,  a  tendency  toward  homogenization  is  anticipated  to  be 
common  in  field  conditions.  Most  wells  should  experience  strong  redistribution  effects,  but 
some  wells  may  maintain  stratification  or  perhaps  re-stratify  differently  from  the  surrounding 
formation.  Variations  in  chemical  concentrations  during  purge  sampling  indicate  that  the  same 
well/aquifer  flow  dynamics  are  at  play  while  purging.  Passive  sampling  does  not  alternative 
mixing  phenomena,  while  purging  tries  to  overcome  these  unknowns  through  well  clearing. 
Both  are  imperfect  in  providing  proof-positive  of  “representativeness”  of  collected  samples. 
However,  this  research  shows  that  under  common  conditions  the  mechanics  of  averaging  by  both 
approaches  is  similar  in  its  end  result.  Ongoing  technical  transfer  of  these  findings  will  promote 
better  understanding  in  the  environmental  community  that  wells  often  represent  a  mixed  flow- 
weighted  average  of  the  adjacent  formation  chemistry.  This  better  understanding  will  yield  cost 
savings  in  both  short-term  and  long-term  timeframes  by  accelerating  the  approval  process  for 
non-purge  alternative  sampling  strategies,  including  passive  sampling  and  in  situ  sensor 
technologies. 
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APPENDIX  A.  SUPPORTING  DATA 


Table  A-1;  2  June  2009  Sampling  Event  Lab  Data,  Zone  Sampling  with  Baffles 

Table  A-2;  17  June  2009  Sampling  Event  Lab  Data,  Zone  Sampling  without  Baffles 

Table  A-3;  24  August  2009  Sampling  Event  Lab  Data,  MW-24  Purge  Event  Complex 

Table  A-4;  25  August  2009  Sampling  Event  Lab  Data,  MW-28  Purge  Event  Complex 

Table  A-5;  1  September  2009  Sampling  Event  Lab  Data,  MW-29  Purge  Event  Complex 

Table  A-6:  2  September  2009  Sampling  Event  Lab  Data,  MW-27  Purge  Event  Complex 

Table  A-7:  16  March  2010  Sampling  Event  Lab  Data,  MW- 15  Purge  Event  Complex 

Table  A-8;  16  March  2010  Sampling  Event  Lab  Data,  MW-27  Purge  Event  Complex 

Table  A-9;  17  March  2010  Sampling  Event  Lab  Data,  MW- 16  Purge  Event  Complex 

Table  A-10;  17  March  2010  Sampling  Event  Lab  Data,  MW-28  Purge  Event  Complex 

Table  A-1 1 ;  16  August  2010  Sampling  Event  Lab  Data,  MW-24  Purge  Event  Complex 

Table  A-12;  20  July  2011  Sampling  Event  Lab  Data,  MW-27  and  MW-28  Tube  Wells 

Table  A-13:  9  August  2011  Sampling  Event  Lab  Data,  MW-27  and  MW-28  Tube  Wells  and 
Baffled  PDB  Sampling 

Table  A-14;  MW-27  Tube  Well/  Baffled  PDB  Data  Reduction  and  Plots 
Table  A-15;  MW-28  Tube  Well/  Baffled  PDB  Data  Reduction  and  Plots 


Table  A-16;  Vial  Pilling  Pour  Tests 
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Appendix  B.  List  of  Scientific/Technical  Publications 


Journal  Publications 

Britt,  Sanford  L.,  Louise  Parker,  James  Martin-Hayden,  and  Mitchell  Plummer,  (in  prep) 
Concentration  Variation  During  Well  Purging  in  Contaminant-Stratified  Aquifers.  Manuscript 
in  development  for  submittal  to  Groundwater  or  Groundwater  Monitoring  and  Remediation, 
2014. 

Britt,  Sanford  L.  and  Louise  Parker,  (in  prep).  The  Effect  of  Changing  Pump  Discharge  Rate 
During  Purging,  Manuscript  in  development  for  submittal  to  Groundwater  or  Groundwater 
Monitoring  and  Remediation. 

Parker,  Louise,  and  Sanford  Britt,  (2012),  The  Effect  of  Bottle  Fill  Rate  and  Pour  Technique  on 
the  Recovery  of  Volatile  Organics.  Ground  Water  Monitoring  and  Remediation,  32  (4);  78-86. 

Martin-Hayden,  James,  Mitchell  Plummer,  and  Sanford  L.  Britt,  (2012),  Controls  of  Wellbore 
Flow  Regimes  on  Pump  Effluent  Composition.  Groundwater,  v  52,  p.  96-104 

Martin-Hayden,  James,  Mitchell  Plummer,  and  Sanford  L.  Britt,  (in  prep).  Characteristics  of 
Thermal  Convection  Within  Groundwater  Monitoring  Wells.  Manuscript  in  development  for 
submittal  to  Groundwater,  2014. 

Conference  Presentations  &  Abstracts 

Britt,  Sanford  L.,  (2009),  Downhole  Sensor  Measurements  Before  and  During  Low  Flow 
Purging:  Dynamics  of  Flow-Weighted  Averaging,  Wednesday  Poster  #207,  SERDP/ESTCP 
Partners  Symposium,  Washington,  DC,  December  1-3,  2009. 

Britt,  Sanford  E.,  and  Louise  Parker,  (2010),  Low  Flow  Purging  and  Sampling  in  a  Stratified 
Aquifer:  Purge  Stability  and  Attainment  of  Flow-Weighted  Averages,  Abstract  A-086, 
Proceedings  of  the  International  Conference  on  Remediation  of  Chlorinated  and  Recalcitrant 

Compounds,  Monterey,  California,  May  24-27,  2012. 

Britt,  Sanford  L.,  James  Martin-Hayden,  Mitchell  Plummer,  Louise  Parker,  Jacob  Gibs,  (2010), 
Well  Purging,  Travel  Times,  and  Discharge  Concentrations:  Where  the  Water  Comes  From  and 
When  it  Arrives,  North  American  Environmental  Field  Conference  and  Exhibition,  Tampa, 
Florida,  January  12-15,  2010. 

Britt,  Sanford  E.,  James  Martin-Hayden,  Mitchell  Plummer,  (2010),  Groundwater  Sampling 
Results  may  Differ  when  Sampling  Methods  Differ:  Methods  and  Mechanisms,  North  American 
Environmental  Field  Conference  and  Exhibition,  San  Diego,  California,  January  10-13,  2010. 
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Britt,  Sanford  L.,  and  Louise  Parker,  (2011),  Testing  VO  A  Vial  Pour  Methods  and  Fill  Rates— 
"Slow  and  Smooth"  Convention  is  All  Wrong,  Tuesday  Poster  #159,  SERDP/ESTCP  Partners 
Symposium,  Washington,  DC,  November  28  -  Deeember  1,  2011. 

Parker,  Louise  and  Sanford  Britt.  (2012),  Effect  of  flow  rate  and  fdl  method  on  the  recovery  of 
VOCs.  North  Ameriean  Environmental  Field  Conference  and  Exposition,  February  7-10  in  San 
Diego,  CA  and  March  13-16  in  Tampa,  EL. 

Britt,  Sanford  L.,  James  Martin-Hayden,  Mitchell  Plummer,  Louise  Parker,  (2012),  Sources  of 
Field  Variability  During  Purging,  North  American  Environmental  Field  Conference  and 
Exposition,  February  7-10  in  San  Diego,  California  and  March  13-16  in  Tampa,  Florida. 

Martin-Hayden,  JM,  MA  Plummer,  SL  Britt,  (2010),  Modeling  Hydrodynamics  of  Well-Bore 
Flow  During  Pumping:  The  Black  Box  of  Groundwater  Sampling  Revealed.  SERDP/ESTCP 
Partners  in  Environmental  Technology  Conference,  Nov. -Dec.  2010,  Washington  DC. 

Martin-Hayden,  JM,  (2011),  The  Effects  of  Ambient  Flow  Regimes  Adjacent  to  and  Within 
Screened  Monitoring  Wells.  North  American  Field  Conference,  San  Diego,  CA,  January,  2011. 

Martin-Hayden,  JM,  (2011),  Revealing  the  Black  Box  of  Groundwater  Sampling:  Effects  of  Bore- 
Hole  Flow  and  Mixing  During  Purging.  University  of  Kansas,  Department  of  Geology 
Colloquium,  May,  2011. 

Martin-Hayden,  JM,  (2011),  Analytical  Modeling  of  Flow  Regimes  Within  and  Surrounding 
Screened  Monitoring  Wells  Under  Ambient  and  Pumped  Conditions.  Geological  Society  of 
America  Annual  Meeting,  Minneapolis,  MN,  October,  2011. 

Martin-Hayden,  JM,  MA  Plummer,  SL  Britt,  (2011),  Flow  Regimes  Within  and  Surrounding 
Screened  Monitoring  Wells  Under  Ambient  and  Pumped  Conditions.  SERDP/ESTCP  Partners  in 
Environmental  Technology  Conference,  Nov. -Dec.  2011,  Washington  DC. 
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Appendix  C.  Laboratory  Reports 

AETL 

Job 


Number 

Date 

pages 

Content 

Lab  Report  Number  1 

52918 

5/28/2009 

5 

EB/TB 

Lab  Report  Number  2 

53071 

6/15/2009 

37 

Multilevel  GW 

Lab  Report  Number  3 

53071a 

8/3/2009 

8 

additional  compounds  reported 

Lab  Report  Number  4 

53235 

7/7/2009 

54 

Multilevel  GW 

Lab  Report  Number  5 

53235a 

7/22/2009 

10 

additional  compounds  reported 

Lab  Report  Number  6 

53389 

7/10/2009 

12 

Soil  for  disposal 
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